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Data and 2D scaling relations for galaxies in Abell 1689: a 
hint of size evolution at z^O.2 
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ABSTRACT 

We present imaging and spectroscopy of Abell 1689 (z=0.183) from the GMOS 
multi-object spectrograph on the Gemini-North telescope and the Advanced Camera 
for Surveys (ACS) on the Hubble Space Telescope (HST). We measure integrated 
photometry from the GMOS g' and r' images (for 531 galaxies) and surface photom- 
etry from the HST F625W image (for 43 galaxies) as well as velocities and velocity 
dispersions from the GMOS spectra (for 71 galaxies). We construct the Kormendy, 
Faber-Jackson and colour-magnitude relations for early-type galaxies in Abell 1689 
using this data and compare them to those of the Coma cluster. We measure the intrin- 
sic scatter of the colour-magnitude relation in Abell 1689 to be CTcmr. =0.054 ± 0.004 
mag which places degenerate constraints on the ratio of the assembly timescale to the 
time available (/3) and the age of the population. Making the assumption that galaxies 
in Abell 1689 will evolve into those of Coma over an interval of 2.26 Gyr breaks this 
degeneracy and hmits /3 to be > 0.6 and the age of the red sequence to be > 5.5 
Gyr (formed at z > 0.55). Without corrections for size evolution but accounting for 
magnitude cuts and selection effects, the Kormendy and Faber-Jackson relations are 
inconsistent and disagree at the 2a level regarding the amount of luminosity evolu- 
tion in the last 2.26 Gyr. However, after correcting for size evolution the Kormendy 
and Faber-Jackson relations show similar changes in luminosity (0.22 ± 0.11 mag) 
that are consistent with the passive evolution of the stellar populations from a single 
burst of star formation 10.2 ± 3.3 Gyr ago (z = 1.8±g°g). Thus the changes in the 
Kormendy, Faber-Jackson and colour-magnitude relations of Abell 1689 relative to 
Coma all agree and suggest old galaxy populations with little or no synchronisation in 
the star formation histories. Furthermore, the weak evidence for size evolution in the 
cluster environment in the last 2.26 Gyr places interesting constraints on the possible 
mechanisms at work, favouring harassment or secular processes over merger scenarios. 
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1 INTRODUCTION 

Scaling relations for early type galaxies (ETGs) are one of 
the simplest tools available to study the formation and evo- 
lution of the present day population. They are particularly 
useful for observing changes in the stellar populations as a 
function of galaxy size, velocity dispersion or mass and allow 
us to test the conventional view that ETGs are composed of 
old populations which f ormed in a single burst of star forma- 
tion some 12 - 1 6 Gyr (lBaaddll958l : iTinslev fc Guimll 19761 : 



iBruzual A.lll983l : iHamiltonlligSSl l.The Faber-Jackson rela 



tion (FJR, iFaber fc Jacksonll 19761 ) relates the central stellar 
velocity dispersion [a) to the absolute magnitude or lumi- 
nosity of a galaxy, 

M = OFJRlogCT + /?FJR. (1) 

Likewise, the Kormendy r elation (KR, iKormendvl 119771 : 
iHamabe fc Kormendvlll987l ) relates the effective radius (-Re) 
to the surface brightness of the galaxy (/Xe), 



(^)e = aanlogRe 4- /3k 



(2) 



Email: rcwh@astro.ox.ac.uk 



The original Kormendy relation was defined between the 
surface brightness at Re and Rs itself. It is now conventional 
to use the average surface brightness within Re against Re- 
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For a de Vaucouleurs light profile (|de Vaucouleurg Il953 ) , 
these are rel ated by a constant factor, but for a Sersic profile 
ijSersidl 19631 ^ this is not the case and the conversion depends 
on the Sersic index n. 

Both the FJR and the KR ar e projections of the 
Fundamental plane of E TGs (FP: iDressler et al] 1 19871 : 
iDiorgovski fc Davig|l987t ). linking 7?e, a and (/ie) (some- 
times expressed in Lqpc~^ as {/)e)- 



\0gRe = applogO- + ^Fp{M)e + 7fp 



(3) 



The FP is tilted compared to the Virial Theorem prediction 
which is thought to be caused by a var i ation in mass-to- 
light (M/L) with mass (JFaber et al.lll987l : lRenzini fc Ciottil 

1 19931 ). 

The colour-magnitude relation (CMR, ISandagelll972l ) 
is a tight relation between the (red) colour of ETGs and 
their magnitude (or luminosity). Defining generic variables 
mi, and rur for the magnitudes in the blue and red filters 
respectively, 

mb—mr= OtcmifMr -f ^CMR. (4) 

where Ai refers to the absolute magnitude. 



1.1 Evolution in slope and intercept of scaling 
relations 

Studying scaling relations at earlier cosmic times allows us 
to study the evolution of the stellar populations. Under the 
assumption that all galaxies evolve uniformly in brightness 
from the changes in stellar populations, one should observe 
a change in the zeropoint, /3. For the FJR, the non-local 
relation becomes 



Mz ^ QFjalogcr^ -I- /3fjr + A/?F,m. 



(5) 



where Mz and a^ are the non-local measurements and 
A/?FjR represents the difference in luminosity from the local 
sample. Similarly for the KR, 



{At)e,z = aiis\ogRe,z -f /?KR + A/3k 



(6) 



Studies of the KR at intermediate redshift show changes 
in A/3kr consistent with ETGs becoming more luminous 
with lookback time and the slope gene ra lly appears 
to be unchanged (IBarrientos et al.l 19961: Pahre et al.l 



La Barbera et all l2003l : iFritz etahl l2005l : iHolden et all 
2005|). Similar conclusions arc drawn from the evolution 



19961: ISchade et al.1 [Tggel. Il997l: IZiegler et al.l Il999l: 



of the FJR (iBender et al.l [1996; Zicglc r et all 12001 



20051: ■ 
Barr et al.l 



12005 



Fritz et ahlbOOSl) a n d the FP (Ivan Dokkum fc Franxlll996al 
van Dokkum et all [l998| : iTreu et al.l I2OOII: iKelson et all 



2OOOI: Ivan Dokkum et ail I2OOII: Ivan Dokkum fc EUisI bool: 



2006 



van der Wei et al 



Fritz et al.l l2005l: iHolden et al.l [20051: Ivan der 'Wei et all 



Moran et al. 20051: J0rgensen et al 



van Dokkum fc van der Marell 



2006 



2007 



20061: van der Marel fc van Dokkum 



2007 : iFritz et al.l [JOtoI : IHolden et all I2OI0I : ISagha et ~ 
20101 ). 

An i mportant use of sca li ng re lations was demon- 
strated by iKodama fc Arimotol ([l997 ) who broke the age- 
metallicity degeneracy (|Worthevlll994l ') by showing that the 



CMR slope is the same in clusters at 2 ~ 0.2-0.4 and must 
originate from a variation in metallicity with luminosity. 



rather than age (in that respect, the cluster ETGs were 
found to be coeval and old). 

Recentl y there has b een deb ate over the obser vation of 
downsizing ([Cowie et al.iil996 : Thomas et al. 20051) and the 
role of selection effects: the study of IHolden et al.l (|2010l ) 
at z ~ 0.8 found no evidence that the tilt of the high- 
z FP chan ges from its loca l value, which imphes tha t the 
findings of fjprgensen et al.1 (|2006l ). iFritz et al.l (|2009l ) and 
ISaglia et al.l 1 20101 ) are biased by selection effects. 

We emphasise the importance of updating local scaling 
relations with the latest technology, techniques and cosmol- 
ogy in relation to observations at higher redshift - the vast 
majority of KR, FJR and FP studies compare high redshift 
observations to t he 15 yr old Gunn r-band observations of 
|j0rgensen et al.l l|l996l ) and obse rvations in different b ands 
are even older (see the review bv lD'Onofrio et al.ll2006l ). Al- 
thou gh potentially a good reference, the SDSS FP parame- 
ters ([Bernardi et al.ll2003h have not been used due to their 
disparity with other work. 



1.2 Scatter in the CMR 



iBower et all (| 19921 . hereafter BLE92) first used the U-V 
Colour-Magnitude Diagram (CMD) of the Coma cluster to 
place constraints on the star formation histories (SFHs) of 
the galaxies. They estimated the scatter in the red sequence 
(RS) of the Coma CMD to be ~0.04 mag, which can then 
be used to infer that either: a) all the galaxies in Coma 
formed at high redshift (z>2) but randomly over the pre- 
ceding interval, or b) the galaxies in Coma formed more 
recently but were highly coeval, thus demonstrating a re- 
markable degree of synchronisation. Given that if all the 
galaxies formed at z~l, the star formation would have been 
apparent with moderately deep surveys and had not been 
observed, the authors concluded that the ETGs in Coma 
formed at z>2 with little synchronisation in the SFHs. 

By consideri ng the effects of mo re complex stellar pop- 
ulations (CSPs). [Bower et al.1 (| 19981 . hereafter BKT98) con- 
firmed the original conclusion of BLE92. They also investi- 
gated the effects of merging on the RS and found that the 
small scatter of the CMR requires that the cluster galaxies 
formed in mass sub-umts not much smaller than half their 
pres ent day mass. 

IStanford et al.1 lll995h. lEllis et al.l (|l997h . 

IStanford et al.1 (|l998l ) and iMei et al.1 (|2009l ). have studied 
many clusters with z < 1.3 to address evolution of the 
CMR scatter. They all found remarkably constant scatter, 
generally consistent with that of Coma with values < 0.1 
mag. This implies that the bulk of the stars in these 
cluster ETGs are very old and formed a t z > 2 with 
littl e synchron isation in the SFHs, although IStanford et al.l 
and Mei et al.l highlight the possibility of progenitor bias. 
iMei et al.l also found evidence that the scatter is less 
for more luminous galaxies, suggesting more massive 
galaxies are sl i ghtly older (by ~ 0.5 Gyr). We note that 
IStanford et al.l (| 19981 ) discuss a gap in the estimates of the 
CMR scatter between 0.1 < z < 0.3 which, together with 
their slightly higher CMR scatter at 2 > 0.3, led them to 
suggest possible evolution of the CMR scatter during this 
epoch. 
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1.3 This study 

The evolution of tiie KR, FJR, FP and CMR appears to 
be consistent witii passive fading of the stellar population 
from a single burst of star formation at z > 1. Similarly, the 
current evidence from analysis of the CMR is also greatly 
in favour of the simple passive fading of ETGs. However, si- 
multaneous study of multiple scaling relations is rarely per- 
formed and it is unclear if they are all quantitatively consis- 
tent. Furthermore, it is important to account for magnitude 
cuts a nd selection effects given the results of iHolden et alj 
l|2010l ). 

In this paper we present data and analysis techniques 
for the cluster Abell 1689. We compare the 2D scaling re- 
lations (CMR, KR and FJR) for AbeU 1689 to those of 
the local cluster Coma, accounting for magnitude cuts and 
selection effects, and specifically testing for consistency: 
if the evolution of different scaling relations is governed 
solely by luminosity evolution of the stellar populations, 
they should all agree. Note that we will present the FP of 
Abell 1689 i n a future paper. Abell 1689 is a massive cluster 
at z=0.183 (jStruble fc Rood 19991) with an X -rav luminos- 



ity of 20.74x10*^^ erg s~^ (JEbeling et al.lll996h making it one 



of the most X-ray luminous galaxy clusters known (the X- 
ray lu minosity for Coma is 7.21x10** erg s~^. lEbeling et al.l 
1 19961 '). The KR, FJR, FP and CMR have not been previously 
studied in regard to the passive fading of ETGs despite it 
being one of the main targets for HST lensing surveys. We 
present colours and magnitudes for 531 galaxies, as well as 
spectroscopic data for 71 galaxies and surface photometry 
for 43 galaxies with which we investigate the CMR, KR and 
FJR. 

Throughout this work, we adopt a WMAP7 Cosmology 
(jKomatsu et aLlbOllH : specifically, we use Hq = 71km s^^ 
Mpc~^, Jim = 0.27 and f^A = 0.73. All quoted uncertainties 
are standard (la, 68%) unless otherwise stated. The struc- 
ture of this paper is as follows. In iJ5] we discuss the data 
reduction of the CMOS imaging, CMOS spectra and HST 
imaging as well as the data analysis techniques for the pho- 
tometry and the kinematics and end with a discussion of 
how we fit the 2D scaling relations. Then in ij3]we present 
our results and discuss the use of stellar population mod- 
els to interpret them in 311 In 33 we present the discussion 
and conclude in ijG] We also present Appendices regarding 
the data reduction and analysis: Appendix A discusses our 
curve-of-growth (COG) technique and the associated error 
analysis; Appendix B discusses how we matched the spectral 
profile of the stellar library to that of the GMOS-N instru- 
ment, and Appendix C presents the techniques used to fit 
and compare the scaling relations of Abell 1689. 



2.1 GEMINI/GMOS imaging 

Details of the seeing limited images of Abell 1689 are given 
in Table [T] The g'-band images have a plate scale of 0'.'1454 
per pixel (2x2 binning) while the r'-band images have a plate 
scale of 0'.'0727 (no binning). The reduction of theGMOS 
imaging was bas ed on that of |j0rgensen et al.l (|2005l ) and 
iBarr et al.l ((20051) but with a few differences. We use custom 
modified versions of the Gemini IRAF data reduction soft- 
ware (vl.l4) and for completeness, we outline our method 
below. 



2.1.1 Preliminary stages: bias subtraction, mosaicing and 
interpolation 

We first processed the GMOS CCD images by performing 
the bias subtraction using the overscan region. A flat field 
image was created using a series of scaled twilight exposures 
(of differing exposure times) which was then applied to the 
Abell 1689 images. The individual chips were mosaiced onto 
a single image using 'nearest' interpolation to avoid corre- 
lating pixels (the error in nearest interpolation is at most 
half the pixel size which in this case is less than 10% of the 
seeing width, see also H2.5|l . All these processing s.ps were 
performed with gireduce, which makes use of gprepare. 



2.1.2 Scattered light removal 

Before mosaicing the images with GMOSAIC and subtract- 
ing the residual scattered light, the background (sky and 
scattered light) levels on each chip were estimated (using 
the modal value) and equalised (individual chip levels sub- 
tracted, mean level added). However, both the g'- and r'- 
band images show a non-uniform scattered light pattern 
across the three CCD chips. 

We produced a sta cked scattered li g ht fra me following 
the recipe described in |j0rgensen et al.l (|20Q5l ). but to pre- 
vent removal of the intra-cluster light (ICL) a 4th order poly- 
nomial surface was fit to the frame, one chip at a time, to 
provide a smooth representatio n of the scattered light w hile 
leaving the ICL present. As in |j0rgensen et al.l (|2005|), we 
scaled the scattered light frame (a polynomial surface in our 
case) by between 0.8 and 1.2 (in 0.05 intervals) and deter- 
mined the best subtraction by eye for each frame. 

The scattered light corrected frames were then aligned, 
scaled and combined using a modified version of the IM- 
COADD routine (allowing multiple pointings to be combined 
into a single mosaic). The default intensity scaling of indi- 
vidual images by imcoadd (given by the header keyword 
RELINT) was found to be unreliable; instead, we fitted a 
Moffat function to a star present in the overlapping point- 
ings of all images, and used the flux within 3 FWHM to 
scale each of the images to a common standard. 



2 DATA REDUCTION AND ANALYSIS 

We now describe the reduction and analysis techniques used 
to measure integrated photometry (from GMOS-N g_G0301 
and r-G0303 imaging, hereafter referred to as g' and r'), 
surface photometry (from HST/ ACS F625W imaging), and 
kinematics (from GMOS-N MOS). 



2.1.3 Flux calibration 



Flux calibr ation of the g'-b and images followed the pre- 
scription in |j0rgensenl (|2009l ) and flnal zeropoint s are li sted 
in Table [H We used st andards from iLandoltl (|l992l) or 
iLandolt fc Uomotd (|2007l ) to estimate the zero points (ZP) 
for each night. 
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Figure 1. A negative image of the galaxy cluster AbcU 1689 observed using the g'- and r'-band filters with the GEMINI/GMOS imager. 
The g'-band image has been degraded to mateh the resolution of the r'-band image. Overlaid is the HST footprint (solid black square) 
and the ROSAT X-ray map (grey, solid lines; contours were shifted by {-0.008,-0.006} degrees in {RA,DEC} to align with the cluster 
core). Coloured squares mark the galaxies in the spectroscopic sample, apart from #119i #793, #1129 which are guiding/alignment 
stars. Note that #160, #457, #892, #932, #983 and #1051 are not in the cluster potential according to their recession velocities (see 
Table [g. 
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Figure 2. The GMOS-N filters transmission curves (g', r') to- 
gether with the HST/ACS F625W transmiss ion curve are shown 
together with a galaxy template spectrum jSantos et al.ll200M 
which has been redshifted to z=0.183. Also s hown for com parison 
is a standard Johnson V-band filter curve l lBessel]||l990l . photo- 
graphic, from) in the rest frame of the galaxy spectrum. 



The g'-data was all taken during photometric condi- 
tions and with good seeing (~ 0'.'6 in the combined image); 
consequently, the images required little relative scaling. The 
photometric standard was taken minutes after the last expo- 
sure, minimising systematic error from atmospheric changes. 
We scaled all images to the last image of the night, which 
was nearest in time to the observation of the photometric 
standard and also had one of the highest relative through- 
puts of all the images (see Table [2]) . 

The r'-band data were taken over two consecutive 
nights, the first of which was not photometric while the 
second was. From the relative throughput of the images, 
it is clear that the flux in the photometric images of the 
second night varies only by a few percent compared to the 
non-photometric images of the first night, and we scaled all 
images to the brightest photometric image of the second 
night (with airmass <1.2) and flux calibrated the images 
using an average of the ZPs derived from all standards. We 
also s caled the science images to correct for airmass (accord- 
ing to I j0rgensenll2OO9t ) . The ZPs in Table [5] do not include a 
term for the extinction from the Galaxy. For t he coordinates 
of Abell 1689, Ac, = 0.10 6 and Ar = 0.075 (|Schlegel et al.l 
119981 : ICardelli et al.lll989l ): these corrections are applied to 
all integrated ( ij2.3p and surface photometry ( i|2.3.ip mea- 
surements. 



2.2 HST imaging 

We used //ST/ ACS images of Abell 1689 obtained with the 
Wide Field Channel (WFC) as part of the ACS Guaran- 
teed Time Observation program 9289 (P.I. H. Ford). The 
WFC detector consists of two 4096 x 2048 SITe CCDs but- 
ted together along their long dimension and separated by a 
gap corresponding to approximately 50 pixels (2.5 arcsec). 
The plate scale is 0.050 arcsec pixel"^. Each chip uses two 
amplifiers to read a single 2048 x 2048 quadrant. We down- 
loaded from the public HST archive the images employing 
the F625W filter, which resembles the SDSS/GMOS r' fiher 



(see Fig. [2]). The observations were carried out on 14 June 
2002 and consist of eight exposures (four pointings split in 
two to deal with cosmic rays) for a total integration time 
of 9500 sec. The dithered images were able to cover the gap 
between the adjacent chips. 

All images were calibrated using the standard reduction 
pipeline PYRAF/CALACS maintained by the Space Tele- 
scope Science Institute. Reduction steps include bias sub- 
traction, dark current subtraction, and fiat fielding, as de- 
scrib ed in detail in the ACS instrument ([Mavbhate A., et al.l 
I2OIOI ) and data (|Pavlovskv C. et al.ll2004l ) handbooks. 

We used PYRAF task MULTIDRIZZLE together with 
standard tasks in IRAFJj to combine all exposures into 
a single geometrically corrected image while rejecting cos- 
mic rays. We let MULTIDRIZZLE generate the inverse- 
variance weight ing map automatically, w hich is suitable for 
SExTRACTOR (JBertin fc Arnoutd 1 19961 ). We analyzed the 
images obtained with different kernels by comparing the 
point spread function (PSF) FWHM derived fitting a two- 
dimensional Gaussian to a few stars distributed in the field 
of view. We chose to use the SQUARE drizzle kernel. In ad- 
dition, for the photometric analysis we need a well sampled 
PSF. The most suitable values for scale and pixfrac were 
found to be 0.030 arcsec pixel"^ and 1.0, respectively. Prior 
to MULTIDRIZZLE, we processed every flat-fielded image 
by fitting a sky level in each quadrant separately, as sug- 
gested by iSirianni et al.l (|2005|), because of residual differ- 
ences after the bias subtraction. 

The final, combined, and geometrically corrected image 
has a size of 7062 x 7251 pixel^, with 0.030 arcsec pixel"^. 
The resulting field of view is approximately 212" x 218" in 
the shape of a rhomboid. 

The transformation to the AB photometric system fol- 
lows [Sinanni_eLalJ ((20051). We adopted the photometric key- 
wo rds tabulated in the header of the images, as suggested 
by IPavlovskv C. et al.l (|2004r ) to ensure up-to-date values 
based on the throughput curves of all HST optical compo- 
nents: 



TAB = -2.5 log/(F625W) -f 25.9186, 



(7) 



where /(F625W) refers to the integrated fiux in units of 
electrons per second. 

We produced a suitable PSF for each sample galaxy in- 
cluded in the HST image. The aim is to obtain realistic PSFs 
that underwent the same reduction-steps we performed on 
the observ ed images. The metho d is similar to the one 
adopted bv lChiboucas et al.l ((20091). For each galaxy we de- 
rived from the flat-fielded images the coordin ates in pixel of 
its ce ntre. By using the TINYTIM package ((Krist fc Hookl 
(1999I) we produced PSFs corresponding to all the exposures. 
We created TINYTIM PSFs with a diameter of 3 arcsec. 
The PSFs were added to synthetic images equal to the flat- 
flelded ones, but with pixel values set to zero. These were 
then combined, aligned, and corrected for geometric distor- 
tion in the same way as the real images. Finally, we extracted 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories which arc operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Filename 



Airmass 

sec(z) 



Moffat 
FWHM 



Gaussian 
FWHM 



Scaling 



Relative 
Throughput 



Filter 



Obs. 
Time (UT) 



Exp. Time 



N20031224S0065 


1.327 


0'.'56 


0'.'619 


1.059 


0.971 


g' 


14:57:29 


180 


N20031224S0066 


1.311 


0'.'53 


0('583 


1.062 


0.964 


g' 


15:01:24 


180 


N20031224S0067 


1.295 


0'.'53 


0^593 


1.060 


0.962 


g' 


15:05:20 


180 


N20031224S0068 


1.281 


0'.'52 


0'.'585 


1.059 


0.960 


g' 


15:09:14 


180 


N20031224S0070 


1.243 


0'.'57 


0'.'597 


1.006 


1.001 


g' 


15:20:09 


180 


N20031224S0071 


1.231 


0"54 


o'.'eio 


1.009 


0.995 


g' 


15:24:04 


180 


N20031224S0072 


1.219 


0'.'52 


o'.'eoo 


1.000 


1.000 


g' 


15:28:00 


180 


N20031224S0073 


1.208 


0'.'50 


0^580 


1.000 


0.998 


g' 


15:31:54 


180 


Coadded g'-band 


1.219 


0'.'56 


0'.'607 


1.000 


1.000 


g' 


- 


1440 


N20011223S174 


1.310 


1'.'13 


1'.'159 


1.005 


1.031 


r 


15:01:09 


300 


N20011223S175 


1.279 


1'.'08 


1'.'164 


1.009 


1.020 


r' 


15:09:12 


300 


N20011223S176 


1.252 


1'.'07 


1'.'147 


1.007 


1.015 


r' 


15:17:14 


300 


N20011223S178 


1.213 


1'.'08 


1'.'118 


1.012 


1.000 


r' 


15:29:56 


300 


N20011223S179 


1.191 


1'.'14 


1'.'246 


1.012 


0.994 


r' 


15:38:10 


300 


N20011223S180 


1.172 


1'.'02 


1'.'066 


1.009 


0.992 


r 


15:46:14 


300 


N20011224S139 


1.518 


1'.'25 


1'.'158 


1.021 


1.071 


r' 


14:17:41 


300 


N20011224S140 


1.466 


1'.'30 


1'.'134 


1.015 


1.063 


r' 


14:25:54 


300 


N20011224S149 


1.167 


0'.'85 


0'.'864 


1.000 


1.000 


r' 


15:44:19 


300 


N20011224S150 


1.150 


0^88 


0'.'997 


1.007 


0.989 


r' 


15:52:22 


300 


N20011224S151 


1.135 


- 


- 


- 


- 


r' 


16:00:28 


300 



Coadded r'-band 



1.167 



1'.'06 



0'.'936 



1.000 



1.000 



3000 



Table 1. The GEMINI/GMOS imaging data. The g'-band data were observed with a plate scale of 0'.'1454 (2x2 binning) while the 
r'-band data were observed with a plate scale of 0^'0727 (1x1 binning). The last exposure in the r'-band was badly affected by twilight 
and was rejected because the scattered light pattern could not be adequately subtracted off. Images arc scaled to the exposure nearest 
in time to the observation of the photometric standard (see Table [2]l 



Filename 



Object 
Name 



Airmass 

sec(z) 



Exp. Time 



Obs. Time 
(HH:MM:SS) 



Counts 

(e-/s) 



mstd 
(mag) 



mzp 

(mag) 



N20031224S0081 
N20031224S0082 
N20031224S0083 
N20031224S0084 



PG1323-086 
PG1323-086 
PG1323-086 
PG1323-086 



1.241 
1.239 
1.237 
1.235 



16:06:09 
16:06:58 
16:07:49 
16:08:36 



6.17x10^ 

6.38x10^ 

6.00e5 

6.24e5 



13.335" 
13.335^,13.283'= 
13.335'", 13.283'= 
13. 335'", 13. 283'= 



27.844 
27.881 
27.814 
27.856 



Filter 



N20011223S119 
N20011224S153 
N20011224S154 



G191B2B 
PG1323-086 
PG1323-086 



1.197 
1.219 
1.216 



09:42:04 
06:09:38 
06:09:38 



24.90x105 
4.93x10^5 
5.25x10^ 



12.044'' 
13. 663'", 13.660'= 
13. 663'", 13. 660'= 



28.062 
27.919 
27.987 



Table 2. The GEMINI/GMOS imaging standards used to calibrate the photometry used here. Note that the g'-band standards were 
observed straight after the science observations, while the r'-ba nd standards wer e observed some hours before hand (see Table [TJ. 
Zero points (ZP) were calculated according to the prescrip tion of [j 0rgcnsc n| ll200gf l. Magnitudes for the photometric standar ds (m gtd) 



were calculated usi ng the { (B, B — V) transformation in ISmith e t al. (20,0^) and the data fro m either Landoltj (suffix a, 1 19921 ) or 
iLandolt &: Uomotd (suffix b. 120071 ): magnitudes for these standards calculated bv |j0rgensenl 1120091 ) are also given (suffix c). The average 
ZP from the above standards for the g'-band data is 27.85±0.03 while the average ZP for the r'-band is 27.99±0.07. 



the PSF from the processed image, with a pixel scale of 0.03 
arcsec pixel~^. 



2.3 Integrated Photometry 

We used SExtractor on the GMOS g'- and r'-band images 
to calculate the integrated magnitudes in those filters. We 
selected galaxies using the g' image (it has better seeing 
than the r' image, see Table [l]) to be above a threshold of 3 
times the sky noise and have a minimum area of 8 pixels. A 
total of 64 de-blending sub-thresholds were permitted with 
a minimum contrast of 1 x 10^^. No filter was applied for 
detection. 



When running SExtractor in dual image mode, the 
images must be identically aligned, scaled, sampled and PSF 
convolved. We re-binned and aligned the r'-band image to 
match the g'-band image and degraded the g'-band resolu- 
tion to that of the r'-band by convolving with a Gaussian 
(FWHM=0'.'87'fl). A pseudo-colour image of Abell 1689 us- 
ing the g'- and r'-band images is presented in Fig. \T\ 

We extract total (MAG_AUTO: based on the Kron def- 
inition) and fixed circular aperture magnitudes (MAG_AP: 
diameter of 2'.'9) for all detected objects. This aperture size 

^ This number is found by trial and error to give similar Gaussian 
fits to both the g'- and r'-band images. 
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was chosen because it is approximately three times the Gaus- 
sian FWHM of the seeing disc in the r' image and exactly 
20 pixels. From the aperture magnitudes in g' and r', we 
calculate the g'-r' colour of the galaxies. 

The errors produced by SExtractor when using just 
the input science images appeared to be underestimated. 
Instead, when provided with a variance image to act as a 
weight map (produced by the GMOS pipeline with stars 
masked), SExtractor produced more realistic errors, used 
here. 

We present the colour-magnitude diagram in Fig. [S] the 
r' magnitude represents the total (Kron) magnitude while 
the g'-r' colour is derived from the aperture magnitudes. 
To use only reliable photometry of galaxies, we select data 
where: the errors on magnitudes are less than 0.1 mag; total 
magnitudes are brighter than 25 mag in g' and r'; SExtrac- 
tor has not flagged the photometry for any reasor|j); the 
object was not identified as a star (class_STAR < 0.99) and 
the object has a FWHM> 4 pixels. Note that the error el- 
lipses neglect the correlation between the r'-band aperture 
and total magnitudes. This data is further presented in H3.ll 
and discussed in i]5.1l We show the constraints for spectro- 
scopic sample selection as dotted lines, which are discussed 
in il31 



In a future paper we will compare these results with full 
surface fitting and decomposition methods. We developed 
our own curve-of-growth (COG) software (written in IDL) 
and implemented this software on the HST/ACS F625W 
image for the spectroscopic sample. This limited the final 
sample which could enter our analysis: we observed a total 
of 77 galaxies in the spectroscopic sample and while all these 
galaxies have ground-based g'- and r'-band imaging ( i)2.ip . 
the HST image covers only 43 of the original 77 galaxies. 

Details of the COG algorithm (including masking of 
nearby sources, correction of PSF effects) are given in Ap- 
pendix|X] We use Monte-Carlo simulations of the fitting pro- 
cedure, to estimate the random and systematic uncertainties 
in logi?e and log/e. These simulations include photon shot 
noise, typical systematic errors in the subtracted sky level 
and a suitable range of radii over which we fit model COGs 
(see Appendix[X|. We calculate the covariance terms for use 
when fitting the scaling relations (see tj2.9p and plotting er- 
ror ellipses in ^ Both the de Vaucouleurs and Sersic COG 
surface photometry are given in Table 21 (/i)e has been cor- 
rected for cosmological dimming and logi?e is given in kpc 
(assuming the cosmology described in iJT]). 



2.3.1 Surface Photometry 

There are two common approaches to calculating the ef- 
fective radius r^ and the average surface brightness within 
re, {^J)e from the fi ux calibrated i mages : curve-of-growth 
fitting (GOG, e.g. [Dressier et all 1 19871 : |j0rgensen et all 
19921) and differential surface fitting (e.g. GASP2D 



Mendez-Abreu et al]|2008l ). The COG technique was origi- 



nally used to parameterise the scale and brightness of galax- 
ies, assuming a de Vaucouleurs like profile for all galaxies, 
where the surface brightness I{R) is defined to be 



I{R) = /oexp ('-7.676[i?/i?e] 



1/4 



(8) 



However, more recently the differential fitting tech- 
nique has been more widely used in the literature (e.g. 
Ivan Dokkum fc Franxl Il996bl ): furthermore it has become 
common to fit the more general Sersic model (for which the 
de Vaucouleurs profile is a special case), defined as 



I{R) = Joexp -b[R/R, 



,1/r 



(9) 



where b ~ 2n — 0.324 over the range 2 < n < 10 
((Ciottil [l99l|) . Recall also that the intensity at Re is Je = 
/o exp(— 6) and the average intensity within Re is {I)e = 
le exp(6)r(2n)n6~^" where F is the gamma function. 

We chose to use the older COG technique for consis- 
tency with the majority of earlier work on scaling relations 
and fit both de Vaucouleurs and Sersic profiles, for which 
the integrated luminosity increases with projected radius as 

L{R) = IoRl^i[b{R/Re)'^", 2n] (10) 

where 7 is the incomplete gamma function (|Ciottil Il991f ) . 



^ other than flags 1 and 2 which inform if the photometry was 
biased by nearby bright objects or originally blended with a neigh- 
bour, respectively. 



2.4 Photometric accuracy 

We perform a number of cross-checks to ensure the accuracy 
of our photometry. The results of these checks are discussed 
below and when (internally) comparing the derived magni- 
tudes of the galaxies we quote the mean difference and the 
RMS scatter in the relation. 

We check the GMOS r'-band ZP in two ways: firstly 
we compare the magnitude of an unsaturated star in our 
field (RA=13:11:33.5, DEC=-01:20:44.7) to the same stars 
observed with SDSS; we measure an r' magnitude of 17.12 
while SDSS quotes an r ' magnitude of 17.150 ± 0.005 (un- 
fortunately the same star is saturated in the g'-band due to 
better seeing). 

We also check how consistent the de Vaucouleurs and 
Sersic COG magnitudes are: a direct comparison in the HST 
F625W image indicated that the Sersic magnitudes were 
0.06 mag fainter than the de Vaucouleurs magnitudes (with 
0.2 mag scatter; the median difference was also 0.06 mag). 

Finally we compare the GMOS r' SExtractor magni- 
tudes to the HST r' COG magnitudes: we find the SExtrac- 
tor Kron magnitudes are fainter than the de Vaucouleurs 
magnitudes by 0.12 mag (with 0.20 mag scatter and median 
of 0.17 mag) and fainter than the Sersic magnitudes by 0.08 
mag (with 0.22 mag scatter and median of 0.13 rnag). We 
estimate that the HST F625W and GMOS r' images are 
calibrated to better than 0.1 mag (the scatter in the r' ze- 
ropoints in Table [2] is 0.07 mag). We believe the remaining 
difference is because SExtractor (Kron) magnitudes are 
not extrapolated to infinity. In principle, we could gener- 
ate simulations to test this, but feel it is not central to our 
analysis because any uncertainty in the GMOS r' calibration 
has little bearing on the results, as we use the scatter in the 
CMR derived from the GMOS r' SExtractor catalogue, 
which is independent of the zeropoint. 
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Table 3. The SExTHACTOR photometry of the GMOS r' and g' images. The Phot. ID is the identification from SExTRACTOR (which 
wc find no need to quote other than in this table), while the Spec. ID is the ID in the spectroscopic sample (which is quoted in this 
paper when referring to a galaxy or star). R.A. and DEC. are given in degrees and are calculated from the g' GMOS image coordinates. 
The Red Seq. column informs us if the object is more likely to be in the red sequence (Y=yes) than not (N=No; i.e. it is more likely 
to bo in the outlier distribution) using the results of the mixture model (see il2.9l ^3.11 & ilC)l . The apparent r' magnitudes and the 
g'-r' colours have been corrected for extinction (atmospheric and Galactic). We refer the reader to the SExTRAOTOR manual for further 
information on the flag codes (v. 2. 13, §9.1, p28). Note that the g'-r' colour is an aperture magnitude while the r' magnitude is a Kron 
(total) magnitude. Note also that galaxy #655 does not meet the criterion specified in i|2.3l but wc include it here for completeness (it 
is included in the KR and FJR). The full table appears in the online version. 



Phot. 


Spec. 


R.A. 


DEC. 


Red 


r' 


g'-r' 


r'flag 


g'flag 


ID 


ID 


(deg) 


(deg) 


Seq. 


(mag) 


(mag) 






154 


14 


197.8058167 


-1.3739161 


Y 


19.2334±0.0014 


0.9509±0.0008 


2 


2 


1364 


15 


197.8058014 


-1.3269255 


Y 


17.7165±0.0007 


0.9600±0.0004 


2 


2 


261 


30 


197.8095551 


-1.3659689 


Y 


18.5313±0.0011 


0.9600±0.0006 


2 


2 


1317 


45 


197.8116913 


-1.3224429 


Y 


20.0542±0.0021 


0.9136±0.0013 








1213 


70 


197.8158112 


-1.3199574 


Y 


18.9672±0.0009 


0.9675±0.0007 


3 


3 


1217 


74 


197.8159485 


-1.3183731 


Y 


18.6034±0.0008 


1.0388±0.0006 


3 


3 


1508 


91 


197.8197632 


-1.3340420 


Y 


19.3286±0.0018 


0.9885±0.0011 


2 


2 


1658 


135 


197.8251495 


-1.3040887 


Y 


19.6441±0.0014 


0.9673±0.0014 








1149 


152 


197.8277740 


-1.3182346 


Y 


19.2382±0.0009 


1.0686±0.0008 








435 


160 


197.8289032 


-1.3585705 


N 


19.5570±0.0016 


0.2652±0.0007 








891 


217 


197.8361206 


-1.3346049 


Y 


18.2879±0.0006 


1.0289±0.0004 


2 


2 


1197 


218 


197.8362885 


-1.3172021 


Y 


20.5914±0.0027 


0.9489±0.0021 








1430 


237 


197.8394165 


-1.3294218 


Y 


18.6605±0.0009 


0.9995±0.0006 


2 


2 


385 


239 


197.8395233 


-1.3621749 


Y 


20.1782±0.0019 


0.9885±0.0014 








471 


259 


197.8419952 


-1.3572158 


N 


19.4536±0.0011 


0.7208±0.0007 








461 


286 


197.8456116 


-1.3554035 


Y 


18.9960±0.0009 


0.9863±0.0007 


2 


2 


1237 


323 


197.8496857 


-1.3192796 


Y 


19.4914±0.0016 


0.9888±0.0011 


2 


2 


520 


341 


197.8513641 


-1.3536860 


Y 


18.6636±0.0007 


0.9992±0.0006 


3 


3 



2.5 GEMINI/GMOS-N spectroscopy 

The GMOS-N instrument in multi-object spectroscopy 
(MOS) mode simultaneously provides multiple spectra in 
the range 0.35^m - 1.1pm inside a 5.5' x 5.5' field of view 
iHooket al.ll2004l ). The spectra of 77 galaxies in Abell 1689 
and 3 guide/alignment stars were taken with GMOS-N be- 
tween 10th January 2002 and 15th January 2002 using 
the B600_G5303 grating and 0'.'75 slits, giving a resolution 
R ~ 1700 (ct ~ 75km s~^) and a (observer-frame) wave- 
length range of approximately 3500A to 7000A, depending 
on the position in the field of view (FOV). The spectro- 
scopic sample is bluer than g' — r' = 1.1 mag, redder than 
g' -r' = -0.17r' -I- 3.4 mag and 50% complete to r' = 19.5 



mag (see gC2|. Of the total number of galaxies matching 
these limits, around 50% are in the spectroscopic sample 
(but note that as a function of magnitude, we are consistent 
with sampling 60% of each magnitude bin, see i]C2[l . The lo- 
cations of the galaxies are shown in Fig.[T]while the selection 
criteria are shown on Fig. [3] The inset histogram in Fig. [T] 
compares the distribution of the parent sample to the total 
spectroscopic sample and the spectroscopic sample that is 
limited to the HST footprint: the total spectroscopic sample 
is a good representation of the parent sample, while the sam- 
ple that lies in the HST footprint is slightly biased towards 
more luminous galaxies. Note that an apparent magnitude 
limit of r' < 19.5 mag is equivalent to a rest frame absolute 
magnitude limit oi Mv < —20.04 mag. 

Four masks were used when observing the cluster and 
for each mask, the data were observed with two different 



central wavelengths (565 nm and 570 nm) so that the two 
chip gaps did not cause discontinuity in the data. The num- 
ber of exposures for each mask varies and details are given 
in Table [S] Our reduction of the spectra is initially very 
similar to that of the imaging ( i]2.ip but differs in e xtrac - 
tion. We again follow the ideas in |j0rgensen et al.l (120051 ) 
and iBarr et al.l (|2005l ) but with a few differences: we use 
custom modified versions of the Gemini IRAF data reduc- 
tion software (vl.l4) and for completeness, we outline our 
method below. 

We process the CCD images in the usual manner, per- 
forming bias subtraction using the overscan region and fiat 
field as normal. The individual chips were mosaiced onto 
a single image using 'nearest' interpolation, which retains 
the statistical independence of the pixels. This has no effect 
on the final result, as GMOS samples the spectral and spa- 
tial PSFs better than the Nyquist limit, by a factor of ~ 5; 
in practice the positioning of our data is accurate to 10% of 
the PSF FWHM. These processing s.ps were performed with 
GSREDUCE, which itself makes use of gprepare, gireduce, 
GMOSAIC and GSCUT. Furthermore, a cosmic ray rejection 
step is present in GSREDUCE. 

Most of the slits were cut at an angle so it was nec- 
essary to rectify the spectra to a uniform space and wave- 
length grid. This was performed with GSWAVELENGTH using 
a 4th-order (Legendre) polynomial fit along A and a 2nd- 
order polynomial fit along the spatial axis. Residuals were 
typically less than 0.2 A. The absolute zero point of the 
wavelength calibration required correction which was ac- 
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Figure 3. The Colour-Magnitude diagram for Abell 1689 created using the GEMINI/GMOS photometry with SExTRACTOR: only data 
with errors less than 0.1 mag are shown. The spectroscopic galaxy sample is highlighted with coloured squares while coloured diamonds 
represent stars. In the legend, labels followed by a star in parenthesis are guiding/alignment stars (#119, #793, #1129); of the other 
diamonds, #972 is reported by SExTRACTOR to be contaminated (the galaxy is peculiar and appears in the g-band image to have two 
parts) while #655 is reported by SExTRACTOR to be a star but its spectrum confirms it to be a galaxy (as does the HST photometry). The 
complete sample is shown by plotting the error ellipses (often too small to distinguish from a dot). Selection limits for the spectroscopic 
sample C il2.5p are shown as dotted lines (r' = 19.5 mag represents 50% completeness) and the inset histogram shows the distribution of 
all objects brighter than r' < 20 for the parent sample (black), for the spectroscopic sample (red) and for the spectroscopic sample inside 
the HST footprint (green). 



complished using the [OI] (5577.34 A) and NaD (5889.95 
A, 5895.92 A) sky lines. Our data were reduced to a log-A 
grid with only one interpolation to help retain the statistical 
independence of the pixels and minimise correlations. 

The sky was removed using GSSKYSUB and a ID galaxy 
spectrum extracted using an aperture of diameter l'.'4 (4.3 
kpc at z=0.183) tracing the peak flux of the galaxy spec- 
trum. Individual ID spectra from each exposure were com- 
bined by scaling by the median, rejecting bad pixels (includ- 
ing the chip gaps) and clipping the data at ±5ct to remove 
bad pixels or cosmics that escaped earlier detection. 



2.6 Kinematics 

To extract the kinematics from the spectra, we assume that 
the galaxy spectrum G (sampled in logA) is the convolution 



of a stellar template T with a distribution of stellar velocities 

L{v), 



G(logA) =T{log\)(g)L{v). 



(11) 



Furthermore, we assume that L{v) takes the parameterised 
form of a Gaussian with velocity V and velocity dispersion 

O". 

In order for Eq.[TT]to be valid, the stellar template pop- 
ulation must match the galaxy population and both should 
have the same spectral resolution (i.e. they should both ap- 
pear to have been observed with the same instrument). In 
practice, this is achieved by taking a large stellar library 
(such as the Indo-U.S. Library of Coude F eed Stellar Spec- 
tra (CFLIB) used here. lValdes et aLlliooi ) observed at high 
spectral resolution and degrading it (via convolution) to 
match the spectral line profile of the instrument; one then 



10 R. C. W. Houghton, Roger L. Davits, E. Dalla Bonta, R. Masters 



Table 4. HST F625W COG surface photometry for 43 galaxies which are also in the GMOS spectroscopic sample. Both de Vaucouleurs 
and Scrsic COGs were used: de Vaucouleurs parameters are headed with (dV) while Sersic parameters are headed (S)\ the Sersic index n 
only appears in this table for Sersic COG fits (n = 4 for de Vaucouleurs COGs). Values of logiJe are in kpc and values of (/i}e are corrected 
for extinction (atmospheric and Galactic), Tolman dimming and the bandpass term of the K-correction (assuming the cosmology described 
in [[TJ. The apparent magnitudes m are calculated directly from the counts in the COG and have only been corrected for extinction, as 
described in i|2.1.3l Uncertainties in the parameters (tr) were calculated from simulation, including the correlation coefficient between 
logRe and (/i)e (p, see ^2.3.1I &: iJAS}. The flags represent the quality of the COG flts; 1 for galaxies whose COG was well represented by 
the model COG; 2 for galaxies whose COG was not well represented by the model COG; 3 for galaxies whose COG did not asymptotically 
tend to a constant value and/or had severe contamination from neighbouring galaxies such that, even after masking, were unreliable (see 
^A2|l . Only galaxies with flags 1 or: 2 enter the fitting process; galaxies #584 and #610 are BCGs which were fitted in some cases (c.f. 
33] and Table 0. 



Spec. 


logiJe 


{t^)e 


P 


ni 


Flag 


logiJe 


(m> = 


P 


m 


n 


Flag 


ID 


(dV) 


(dV) 


(dV) 


(dV) 


(dV) 


(S) 


(S) 


(S) 


(S) 


(S) 


(S) 


286 


0.40±0.35 


19.79±1.10 


-0.998 


18.78±0.00 


2 


0.35±0.13 


19.66±0.39 


-1.000 


18.86±0.00 


2.8 


2 


341 


0.37±0.35 


19.34±1.10 


-0.999 


18.47±0.00 


2 


0.35±0.13 


19.28±0.39 


-1.000 


18.52±0.00 


3.2 


2 


368 


1.02±0.34 


21.37±1.10 


-0.999 


17.24±0.00 


2 


0.85±0.03 


20.77±0.07 


-0.999 


17.50±0.00 


2.0 


2 


371 


0.40±0.35 


20.38±1.10 


-0.998 


19.33±0.00 


1 


0.47±0.20 


20.65±0.71 


-1.000 


19.24±0.00 


5.2 


1 


390 


0.41±0.35 


19.74±1.10 


-0.998 


18.65±0.00 


1 


0.37±0.13 


19.59±0.39 


-1.000 


18.71±0.00 


3.2 


1 


398 


0.25±0.35 


18.87±1.10 


-0.998 


18.60±0.00 


2 


0.22±0.11 


18.80±0.29 


-1.000 


18.67±0.00 


2.6 


1 


433 


0.70±0.37 


22.11±1.13 


-0.993 


19.57±0.00 


2 


0.57±0.20 


21.68±0.46 


-1.000 


19.77±0.00 


2.0 


2 


435 


0.63±0.35 


20.62±1.10 


-0.999 


18.42±0.00 


1 


0.68±0.08 


20.79±0.29 


-1.000 


18.36±0.00 


4.9 


1 


463 


0.60±0.35 


20.45±1.10 


-0.999 


18.43±0.00 


1 


0.63±0.08 


20.57±0.29 


-1.000 


18.38±0.00 


4.8 


2 


476 


0.44±0.35 


19.17±1.10 


-0.999 


17.95±0.00 


2 


0.39±0.05 


19.03±0.15 


-1.000 


18.03±0.00 


2.8 


2 


481 


0.25±0.35 


19.31±1.10 


-0.998 


19.02±0.00 


2 


0.26±0.11 


19.40±0.29 


-1.000 


19.08±0.00 


2.4 


1 


501 


0.82±0.35 


21.01±1.10 


-0.999 


17.86±0.00 


2 


0.80±0.06 


20.92±0.17 


-1.000 


17.91±0.00 


3.3 


2 


508 


0.57±0.35 


20.16±1.10 


-0.999 


18.30±0.00 


2 


0.43±0.06 


19.70±0.14 


-0.999 


18.50±0.00 


2.2 


2 


514 


0.45±0.35 


19.83±1.10 


-0.998 


18.57±0.00 


2 


0.94±0.09 


21.80±0.33 


-1.000 


18.06±0.00 


9.8 


2 


531 


0.35±0.35 


20.27±1.10 


-0.998 


19.48±0.00 


2 


0.33±0.20 


20.29±0.46 


-1.000 


19.59±0.00 


1.8 


2 


549 


0.66±0.35 


20.45±1.10 


-0.999 


18.09±0.00 


1 


0.64±0.06 


20.35±0.17 


-1.000 


18.13±0.00 


3.5 


1 


567 


0.12±0.35 


18.64±1.10 


-0.998 


18.99±0.00 


1 


0.13±0.16 


18.67±0.51 


-1.000 


19.00±0.00 


3.6 


1 


584 


1.48±0.34 


22.07±1.09 


-0.999 


15.62±0.00 


3 


0.99±0.00 


20.47±0.01 


-0.998 


16.48±0.00 


1.2 


3 


593 


0.25±0.35 


19.55±1.10 


-0.998 


19.27±0.00 


1 


0.25±0.16 


19.58±0.51 


-1.000 


19.27±0.00 


3.7 


2 


601 


0.56±0.35 


21.01±1.10 


-0.998 


19.16±0.00 


2 


0.46±0.13 


20.65±0.39 


-1.000 


19.30±0.00 


2.9 


1 


610 


1.04±0.34 


20.99±1.10 


-0.999 


16.77±0.00 


3 


1.05±0.03 


21.03±0.08 


-1.000 


16.76±0.00 


3.9 


3 


635 


0.63±0.35 


20.25±1.10 


-0.999 


18.05±0.00 


3 


0.78±0.09 


20.83±0.33 


-1.000 


17.87±0.00 


5.9 


3 


636 


0.59±0.35 


20.68±1.10 


-0.998 


18.72±0.00 


3 


0.70±0.31 


21.13±1.14 


-1.000 


18.61±0.00 


6.2 


3 


645 


0.60±0.35 


20.02±1.10 


-0.999 


17.97±0.00 


2 


0.55±0.04 


19.83±0.11 


-1.000 


18.06±0.00 


2.6 


1 


655 


0.67±0.35 


20.75±1.10 


-0.999 


18.36±0.00 


2 


0.72±0.07 


20.94±0.23 


-1.000 


18.30±0.00 


4.7 


2 


670 


0.39±0.35 


19.36±1.10 


-0.999 


18.37±0.00 


1 


0.36±0.05 


19.24±0.15 


-1.000 


18.43±0.00 


3.2 


1 


677 


0.81±0.35 


20.98±1.10 


-0.999 


17.89±0.00 


1 


0.92±0.08 


21.38±0.29 


-1.000 


17.76±0.00 


5.2 


1 


690 


0.76±0.34 


19.88±1.10 


-0.999 


17.07±0.00 


2 


0.67±0.01 


19.60±0.03 


-0.999 


17.23±0.00 


2.0 


1 


698 


0.43±0.35 


20.34±1.10 


-0.998 


19.19±0.00 


2 


0.37±0.06 


20.19±0.14 


-0.999 


19.31±0.00 


2.0 


2 


717 


0.33±0.35 


19.29±1.10 


-0.998 


18.61±0.00 


2 


0.33±0.13 


19.33±0.39 


-1.000 


18.64±0.00 


3.1 


1 


723 


0.51±0.35 


19.95±1.10 


-0.999 


18.38±0.00 


2 


0.48±0.05 


19.85±0.15 


-1.000 


18.44±0.00 


3.0 


1 


724 


0.44±0.35 


20.41±1.10 


-0.998 


19.17±0.00 


3 


0.51±0.12 


20.67±0.39 


-1.000 


19.10±0.00 


4.1 


3 


726 


0.63±0.35 


19.90±1.10 


-0.999 


17.73±0.00 


1 


0.65±0.06 


20.01±0.20 


-1.000 


17.70±0.00 


4.0 


2 


753 


0.83±0.34 


19.85±1.10 


-0.999 


16.68±0.00 


1 


0.81±0.03 


19.79±0.08 


-1.000 


16.70±0.00 


3.8 


1 


755 


0.33±0.37 


20.86±1.13 


-0.993 


20.16±0.00 


2 


0.29±0.35 


20.72±1.02 


-1.000 


20.24±0.00 


2.8 


1 


756 


0.95±0.34 


20.61±1.10 


-0.999 


16.81±0.00 


2 


1.19±0.05 


21.53±0.18 


-1.000 


16.53±0.00 


6.5 


2 


814 


0.72±0.35 


20.52±1.10 


-0.999 


17.89±0.00 


2 


0.58±0.04 


20.04±0.11 


-1.000 


18.09±0.00 


2.3 


2 


816 


0.29±0.35 


19.75±1.10 


-0.998 


19.27±0.00 


1 


0.27±0.13 


19.69±0.39 


-1.000 


19.33±0.00 


2.9 


1 


848 


0.41±0.35 


19.64±1.10 


-0.998 


18.55±0.00 


1 


0.40±0.16 


19.59±0.51 


-1.000 


18.58±0.00 


3.4 


1 


852 


0.19±0.35 


19.45±1.10 


-0.998 


19.46±0.00 


2 


0.18±0.22 


19.43±0.57 


-1.000 


19.52±0.00 


2.6 


2 


874 


0.60±0.35 


19.94±1.10 


-0.999 


17.92±0.00 


1 


0.56±0.05 


19.80±0.15 


-1.000 


17.98±0.00 


3.2 


1 


883 


0.81±0.35 


21.18±1.10 


-0.999 


18.10it0.00 


2 


0.63±0.03 


20.59±0.07 


-0.999 


18.42±0.00 


1.2 


2 


906 


0.46±0.35 


20.71±1.10 


-0.998 


19.39±0.00 


2 


0.39±0.20 


20.55±0.46 


-1.000 


19.55±0.00 


1.5 


2 



generates a best-fit template in parallel to fitting the kine- 
matics. 

Usually, both the stellar library and the galaxy spectra 
have spectral profiles closely approximating a Gaussian and 
so convolving with another Gaussian completes the process 
of matching the spectral resolutions. However, in the case of 



GMOS the spectral profile is not Gaussian; thus matching 
the spectral profiles was somewhat involved and we dedicate 
an Appendix ([B|) to explaining our technique. 

To calculate the recession velocities V and velocity dis- 
persions a of the individual galaxies, we made use of the 
freely available PPXF (penalised pixel fitting) IDL soft- 
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Mask Exp. Times Nuiii of Frames. Total exposure time Dates of exposures 



Notes 



03 2400 

04 2400,2580 

05 2400,3000 

06 2400,2100,3000 



14400 10/01/2002, 11/01/2002, Light clouds on two exposures taken on 10th 

10800 13/01/2002 

9960 12/01/2002 

15000 14/01/2002, 15/01/2002 Guiding lost on last lOmins of one exposure 



Table 5. Details of the GMOS/Gemini ground based spectra of the galaxies in Abell 1689. 
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are unreliable and we replace them with an upper limit of 50 
km s~^ and do not use them to fit the scaling relations (see 
!j3]). This affects only 4/77 of galaxies, two of which could 
enter the KR and FJR as they have good HST surface pho- 
tometry (#883 & #906). 

We corrected our velocity dispersion measurements 
for aperture effects. Each spectrum was extracted from a 
rectangular aperture of size l'.'4 by I'.'O and corrected to 
the equivalent c ircular radius following the procedure in 
|j0rgensen et al.l (|l995bl ). We also correct our aperture dis- 
persion measurements dap to a velocity dispersion measured 
within a standard apertu re size (1.62 kpc, equiv alent 3''4 
at the di stance of Comaj |j0rgensen et al.lll995a ). The re- 
lation in |j0rgensen et al. (Il995br) is sl ightly different to the 
one found bv ICappellari et al.l ( 20061 ): however, the effect 
of this difference o n our aperture c orrections is small and 
we persist with the |j0rgensen et al.l relation for comparison 
with the literature. 
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Figure 4. A comparison of the velocity dispersions a derived 
from different wavelength ranges: 3500 < A(A) < 4050, 4050 < 
A(A) < 4650, 4650 < A(A) < 5650 and 5650 < A(A) < 7000. 



ware (jCappellari fc Emsellemll2004h to fit the parameterised 
stellar kinematics, and the GA NDALF (Gas AN D Absorp- 
tion Line Fitting) IDL software (jSarzi et al.ll2006r ) to fit any 
emission lines. We divided each input spectrum in four parts 
and fit them separately: wavelengths less than 4050 A; be- 
tween 4050 A and 4650 A; between 4650 A and 5650 A, and 
wavelengths greater than 5650 A (all wavelengths quoted in 
the rest frame of the cluster). We also excluded regions af- 
fected by telluric absorption from the fit, such as the range 
6850 < A(A) < 6950. For each wavelength range a new best- 
fit stellar template was found. Fig. |4] shows the dispersions 
calculated from the middle two sections (4050 A - 4650 A 
and 4650 A - 5650 A) have the least scatter; thus we aver- 
age the dispersions from these two regions for all subsequent 
analysis. 

Although PPXF can determine dispersions less than 
the instrumental resolution (« 75 km s~^ at 0.5 ^m), the er- 
rors (both random and systematic) increase. Consequently, 
we decided that dispersions found to be less than 50 km s^^ 



2.7 The local comparison sample 

To draw conclusions about the changes in the properties of 
the galaxies in Abell 1689 (at z=0.183), we need to compare 
them with local {z ~ 0) galaxies. For a fair comparison, we 
must select the local galaxies from a similar environment 
to the Abell 1689 galaxies (i.e. from a massive cluster) and 
the only suitably dense nearby envir onment is the Coma 
M ouldl 19921) at a distance of 103± 



cluster ( z = 0.024, |H an 



lO Mpc (iLiu fc Grah am 2001, iThomsen et all 1 19971 but see 
also Ijensen et al.l Il999l ) . 

Fortunately, the ETGs in the Coma cluster have 
been rigorously investigated in the past. Like previous au- 



thors (|Ziegler et al.ll200ll : lLa Barbera et al.ll2003l:lBarr et al.l 
[2OO5I), we start with the data from I j0rgensen et al.l (|l999al . 
hereafter J99): a total of 116 velocity dispersions for galaxies 
in Coma were compiled from the literature and homogenised 
in thi s work. We then compare the sample with the photome- 
try of |j0rgensen et al.l (|l995al . hereafter J95a) and find that 
all but 2 galaxies have Gunn r surface photometry avail- 
able (calculated by fitting de Vaucouleurs curve-of-growth 
models). We adopt the formal uncertainties given in J95a 
and J99 but refrain from imposing covariances. Further- 
more, we do not attempt to update the Galactic extinction, 
K-corrections or cosmological corrections but use the data 
as presented and model them (Q as though they were ob- 
served at z = 0.0. We determine that this sample is 50% 
complete to r' = 14.8 (equivalent to M^i = —20.28 ,see 
El. 
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Table 6. The kinematics of the spectroscopic sample within 
±6000km s~^ of z=0.183; out of 80 objects targeted, 3 were guide 
stars leaving 77 galaxies of which 71 were found to meet the cri- 
terion on Kec to be considered part of the cluster. Those galaxies 
not in the cluster (and not shown here) were #160, #457, #892, 
#932, #983 and #1051; of these only #932 has V.ec relative to 
z=0.183 which could be in the cluster (8822km s~^). However, 
none of these galaxies are found to be in the RS according to the 
mixture model (see Table [Sjl. Velocity dispersions measured to be 
< 50km s~^ were deemed unreliable (well below the resolution of 
the instrument) and replaced with an upper limit of 50km s~^; 
these points are not used when fitting the scaling relations, but 
are shown (as diamonds) in Fig. [6] Velocity dispersions given here 
have not been aperture corrected (sec i|2.6|l . The full table appears 
in the online version. 



Spec. 


Kcc 


O-ap 


ID 


(kms-^1) 


(kms-i) 


14 


1672.2±2.4 


83.2±3.5 


15 


-981.8±3.3 


198.4±3.7 


30 


1638.8±2.0 


125.0±2.5 


45 


300.8±4.0 


104.3±5.4 


70 


1439.5±2.7 


134.8±3.3 


74 


2524.8±3.5 


159.1±4.2 


91 


2897.4±5.2 


114.8±6.6 


135 


948.2±2.8 


85.8±4.2 


152 


1291.9±2.9 


158.2±3.4 


217 


299.5±2.9 


172.6±3.3 


218 


2266.3±9.3 


108.6±12.3 


237 


-934.2±4.2 


170.8±4.8 


239 


1545.2±6.3 


146.8±7.6 


259 


1284.4±2.3 


73.6±3.7 


286 


2255.1±2.5 


150.7±2.9 


323 


-83.9±2.9 


90.3±4.4 


341 


802.2±5.1 


201.6±5.5 


343 


-1257.0±4.3 


215.3±4.7 



2.7.1 Companng Coma and Abell 1689 

We wish to infer ages and formation timescales of the 
ETGs in Abell 1689 and Coma by assuming that galaxies 
in Abell 1689 will evolve into those in Coma over 2.26 Gyr. 
This assumption is not immediately obvious; furthermore, if 
untrue it could invalidate our analysis. In a hierarchical pic- 
ture, one expects more massive DM halos to collapse earlier 
and potentially produce galaxies earlier; thus, if Abell 1689 
were considerably more massive than Coma was at z=0.183, 
it is not unreasonable to claim that the Abell 1689 pop- 
ulation would be older, leading to a smaller CMR scatter 
and A/3, thus biasing our measurements. We now look at 
these two clusters in detail and discuss how similar or differ- 
ent they are. We address the X-ray luminosities, the cluster 
masses and the known details of the galaxy populations to 
gauge how robust our approach is. 

While Coma is relatively X-ray bright (Lx =7.21 xlO'''* 
erg s~^), its luminosity somewhat dwarfed by the 
(exceptionally) high X-ray luminosity of Abell 1689 
{Lx =20.74x10** erg s~^). On X-ray luminosity evidence 
alone, one might conclude that the mass of Abell 1689 
is many times the mass of Coma, or that Abell 1689 
has recently undergone a merger (it hasn't, according to 
the relaxed nature of the X-ray contours, iLemze et al.l 



[2003). However, the X-ray luminosity is only part of the 
picture; indeed, detailed X-ray modelling of Coma and 
Abell 1689 indicate that their total masses are quite com- 
para ble given the modelling uncertainties (around 5 x 10^* 



Mq: iMason fc Myers! I2OO0I: lEttori et all I2OO2I: IPeng et all 



I2OO9I ; lRiemer-S0rensen et al.l |2009| ; iMahdavi et al.l l20(m 
The lensing view on the two cluster masses is broadly sim- 
ilar: although Abell 1689 appears marginally more massive 
than Coma, there is considerable scatter in the lensing esti- 
mates (mass estimates for both are around 1 — 2 x 10^ ^ M 



Kubo et al. 200 



Broadhurst et 
2008h . 



)0^; iGavazzi o t al.1 |2009| ; 'Okabc et all I2OI0I; 
ay [2005; Hal kola et al.l 12006: M ahdavi et all 



We note (like many others) that the X-ray masses of the 
inner regions of clusters are genera lly estimated to be half 
those of the le n sing m easurements. lOguri et al.l (|2005f ) and 
iMorandi et al.l (|201ll ) were able to resolve this discrepancy 
with triaxial halos and suggest Abell 1689 is viewed along 
its major axis. This is unlikely for Coma as we see the two 
central BCGs well separated on the sky. Thus certain dis- 
crepancies between Coma and Abell 1689 could in part be 
down to projection effects. 

In terms of mass evolution, it is unlikely that Abell 1689 
will become much more m assive between z=0.183 and now: 
according to the work of iFakhouri et al.l (|2010l ). for halo 
masses > lO^^Mo, the rate of change of mass at < « < 0.2 
is < 4x 10*Moyr~^: therefore, the expected increase in mass 
over 2.26 Gyr is < 9.2 x IO^^Mq, which is a small fraction 
of the total mass. 

Thus we conclude that, aside from projection effects 
from possible triaxiallity, the two clusters are well matched 
in mass and are both reasonably relaxed. 



2.8 Magnitude conversions, Cosmological 

corrections and evolution measurements 

We cannot compare the Coma data to the Abell 1689 data 
directly: the data are in different filter bands using differ- 
ent magnitude systems. Furthermore, the Abell 1689 data 
are affected by cosmological dimming (the Coma data has 
been corrected for this, see H2.7|l and evolution of the stellar 
population (the magnitude of which we wish to determine). 
Therefore, we must address these differences and where ap- 
propriate apply corrections or model the effects (see Q. 

2.8.1 Magnitude systems 

The Coma data of JF95a was observed in the Gunn r-band 
using the Gunn photometric system (relative to a subdwarf 
F6 star) whereas the HST F625W observations are in the 
AB system. We convert the Gunn r photometry to the A B 
system using the corrections listed in lFrei fc GunnI ([199J). 



2.8.2 Cosmological Corrections 

When we observe galaxies at significant redshifts, it is nec- 
essary to correct for cosmological effects (i.e. e xpansion and 
reds hift). We adopt the approaches described in lHogd (|l999l ) 
and iHogg et al.l (120021) to correct our measurements. We 
split the K-correction (JHogg et al.ll2002h into two terms. 



K^Kb + K^ 



(12) 



Data and 2D scaling relations for galaxies in Abell 1689: a hint of size evolution at z^O.2 13 



The bandpass term K}, is easily corrected in the AB magni- 
tude system by reducing the observed brightness by {1 + z). 
However, the colour term Kc depends on the details of the 
underlying stellar population which change with age, metal- 
licity, initial mass function (IMF), dust content, etc.. J95a 
Eissumed no evolution of the galaxy SED when calculating 
the colour term of the K-correction Kc for the Coma data, 
which is a good approximation at low redshift. But at higher 
redshifts, we do not know the exact SED so we do not ap- 
ply a Kc correction to our measurements of Abell 1689. In 
^ we describe how we model our observations by calcu- 
lating magnitudes based on stellar population models in 
blueshifted filter profiles, thus negating the need for Kc- 
However , it is necessary to correct fo r Tolman dimming 
l|Tolman|[l930l : iLubin fc Sandagell200lh . Incorporating this 
and Kh, the bandpass limited AB surface brightness /Xobs 
is related to the bandpass limited rest-frame AB surface 
brightness ^,cst by 

ProBt = A*obB + 7.51og(l + z). (13) 

Furthermore, we can calculate apparent magnitudes using 

m = (M)c,,cst - 51og(i?e) - 2.51og(27r) + 7.51og(l + z) (14) 

and (bandpass corrected) absolute magnitudes using 

M = m-F2.51og(l-h2)-51og(Dz,)-H5 (15) 

where Re is measured in arcseconds and Dl is the luminosity 
distance in pc. 

2.8.3 Size evolution 

Recent literature suggests that galaxies ( both disks and 
spheroids) were more compact in the past (Mo et al.l 19981: 
Ferguson et all l20o4 iBouwens et all l2004l: iTruiiUo et all 



20071 : Ivan der Wei et al.ll2008l : lHopkins et al.|[2009D . The de- 
tails are still very much a topic of debate: the degree of 
comijactncss could depend on galaxy mass |Barden et al.l 
[2005: Mcintosh ct al. 200^; iTruiillo et all 120071 ). We mea- 
sure Re,{fJ,}e and a (for Coma and Abell 1689) with and 
without correcting size evolution. The corrections we ap- 
ply assume that the effective radius now. Re, is related to 
the effective radius at non-zero redshift, Re,z, by a simple 
power-law scaling, 

I0g(i?e)=l0g(i?e,z)+Cl0g(l + Z) (16) 

and therefore, assuming no other changes, 

{M)e = {M)e,.+5Cl0g(l-f2). (17) 

The value of C is still debated: for massive, high concentra- 
tion galaxies. iBouwens et al.l (|2004 ) find C = 1.05±0.21 over 
the ra nge 2.5 < 2: < 6 which agrees with Ivan der Wei et al.l 
(|2008h who found ( = 0.98±0.11 betwe en < z < 1 for mor - 
phologically selected ETGs. However, ITruiillo et"all (120071 ) 
found slightly stronger evolution for massive spheroid-like 
(highly concentrated) galaxies: a factor of 4 ± 0.4 since 
z = 1.5, equivalent to ^ = 1.6. We adopt a value 1^ = 1.0 in 
this work. 

We also calculate the equivalent velocity dispersion 
now, a, from the velocity dispersion measured at higher red- 
shift, az, to be 



with rj — 1.0. In dynamical models of Sersic-like 
galaxies, the projected line-of-sight velocity dispersion 
changes in this w ay whe n Re is scaled, if the Sersic in- 
dex is unchanged ICiottil (|l99ll ). However, iHopkins et al.l 
(|2010|) suggest the Sersic index does change. The effect 
of size evolution on the velocity disp ersion is difBcult 
to me asure and not w e ll con st rained ([van der Wei et al.l 



20091: Ivan Dokkum et all |2009| : 



20081: iCappellari et all _. .___ 

Hopkins et al.1 l2010t ). ICenarro fc TruiiUd (|2009l ) report 



log(o-) = log{a:,) - O.5?7log(l -j- z). 



(18) 



value of 7; ~ .6 w hile the single observation by 
Ivan Dokkum et al.l ((20091) would suggest 77 « 2. In the ab- 
sence of a clear measurement of the change in the Sersic 
index, we adopt -q = 1.0. 



2.9 Fitting 2D scaling relations 

There are many different techniques available to fit (thereby 
represent) scaling relations. The disadvantages and limita- 
tions of common techniques are discussed by iHogg et al.l 
(|2010h who advocate Markov Chain Monte Carlo (MCMC) 
methods. We follow these guidelines and use MCMC meth- 
ods to explore the posterior distributions of the parameters 
in each assumed model; details are given in Appendix [Cl 
We investigate three different models: one where the data 
(with uncertainties) are drawn from a linear relation with 
unknown slope, intercept and intrinsic dispersion and we do 
not attempt to correct for magnitude cuts or selection effects 
(hereafter the linear model or LM); a double linear model 
where two sets of data (with uncertainties) are assumed to 
have the same slope and intrinsic dispersion but with differ- 
ent intercepts and we apply corrections for magnitude cuts 
and selection effects (hereafter the double linear model or 
DLM); and finally, a model where the data are drawn from 
a mixture of two distributions: one being the linear model 
as before, and another outlier distribution having unknown 
first and second moments (hereafter the mixture model, or 
MM). In all cases, we assume normal distributions for the 
models. 

We use the LM and DLM to fit the KRs and FJRs of 
Coma and Abell 1689 to measure any offset (i.e. evolution) 
between them. We use the MM to isolate and measure the 
intrinsic scatter of the CMR in Abell 1689. 



3 RESULTS AND ANALYSIS 

We present the Faber- Jackson relation, the Kormendy re- 
lation and the colour magnitude diagram for Abell 1689, 
below. As discussed in i\1.3\ we will present and analyse the 
FP in a future paper. 



3.1 g'- r' Colour Magnitude Diagram 

The raw g'-r' colour magnitude diagram (CMD) for 
Abell 1689 was presented earlier in Fig. (3] showing only 
the data where the uncertainties were less than 0.1 mag. 
The RS is clearly visible, along with a blue cloud. We also 
over plot (dotted lines) the selection limits for the GMOS-N 
spectroscopic sample and show as an insert the luminos- 
ity distribution of the parent sample and the spectroscopic 
sample. 
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Figure 5. CMD with RS shown (red) and outliers (green) re- 
jected by the MCMC mixture model. 



cessively red objects are found beyond r' > 20 mag, so they 
do not enter our spectroscopic sample. Galaxies redder than 
the RS could be excessively old, excessively metal rich, heav- 
ily extincted and/or intermediate or high redshift interlop- 
ers. Alternatively, we could be underestimating the scatter 
in the CMR at this magnitude (particularly because we as- 
sume a single scatter across all magnitudes'). iTerlevich et al.l 
(|l999l ) correlated the residuals of the CMR with spectral ab- 
sorption line indices and showed that the galaxies scattered 
blueward of the mean relation have an increased hydrogen 
Balmer absorption and are thus younger: if the reverse is 
true, galaxies redder than the RS could be old relics from 
the initial cluster population. It is also worth noting that 
for ages approaching 14Gyr and metallicities around twice 
solar, both the BC03 and M05 models predict a maximum 
g'-r' colour of < 1.4. Thus the few objects redder than g'-r' 
= 1.4 are unlikely to be part of the cluster unless they are 
extremely metal rich or extincted; conversely, the majority 
of the red objects that have g'-r' < 1.4 but are redder than 
the RS could be cluster members, in principle. 



3.1.1 Measured scatter of the CMR in Abell 1689 

We use a mixture model (see i 32.9l fc Appendix[C| to measure 
the slope, intercept and intrinsic scatter of the CMR (Eq. 
|4]). We restrict the galaxies in the sample to be brighter 
than r' < 22 mag as well as the constraints described i]2.3l 
We do not apply constraints on the recession velocities or 
morphologies of the galaxies, which can only increase the 
measured scatter. 

The results of fitting the MCMC mixture model are 
given in Table [TJ in particular, the intrinsic scatter is found 
to be o-cMR =0.054 ± 0.004 mag. 

In Fig. [5] we show the CMD of AbeU 1689 with the 
galaxies marked according to which of the two distributions 
they were most likely to be members of. We also show the 
best fit to the CMR, marginalised over all other parame- 
ters. There is an excess of galaxies below the RS, and a few 
above, which are considered to be outliers by the mixture 
model (points highlighted in green), but the mixture model 
has successfully isolated the RS (points highlighted in red) 
without any cut in colour or clipping imposed a priori. The 
limit in magnitude (r' < 22 mag) affects the resulting scat- 
ter significantly: brighter (lower) magnitude cuts decrease 
the scatter while fainter (higher) magnitude cuts increase 
the scatter. Magnitudes are taken from the entire CMOS 
r-band image, with dimensions 280" x 300". 

BKT98 quote Coma's CMR scatter inside R < 600 kpc 
to be 0.049 for all galaxies brighter than Mv < —18.2 mag 
(98 galaxies). We fit the CMR to galaxies brighter than 
Mv < —17.90 mag inside R <570 kpc. The selection effects 
and resulting scatter in both clusters is very similar; the 
scatter in Coma is slightly smaller t hough uncertainties are 
not given by BKT98. More recently, ITerlevich et all (|200ll ) 
performed a more extensive study of the CMR in Coma: for 
all galaxies brighter than Mv < —18.2 mag in roughly a 
square degree (175 galaxies), they quote an intrinsic scatter 
of 0.063±o;oi* mag, which is statistically consistent with our 
result for Abell 1689 (although our uncertainties are a factor 
of 2 smaller). 

We highlight a handful of objects at relatively faint 
magnitudes (r' ~ 21 mag) that lie above the CMR. These ex- 



3.2 The Faber-Jackson relation 

Table [7| presents the results for fits to the FJRs (Eq. [T]) of 
Coma and Abell 1689 using individual linear models (LMs) 
and the double linear model (DLM). 

In Abell 1689, the calculation of the surface photometry 
for the central BCGs (brightest cluster or cD galaxies: #584 
& #610) was compromised by contamination from other 
galaxies and the ICL. Similar arguments could be made for 
the two BCGs in Coma (GMP2921 & GMP3329). Thus we 
fit the data both with and without the two BCGs and with 
and without correcting for size evolution (see i]2.8.3[) . 

The values for A/3f.tr in Table [7] were either calculated 
by subtracting the absolute magnitude of the F JR at logcr = 
2.2 (where the FJRs are best constrained in individual LMs) 
or directly when fitting the DLM to both the Abell 1689 and 
Coma data. For nearly all cases, the galaxies in Abell 1689 
are fainter in the F625W band than the Coma galaxies in the 
r'-band. The evolution measured from the LMs and DLM 
generally agree, suggesting little bias from cuts or selection 
effects in the FJR. 

Correcting for size evolution has a significant effect in 
decreasing the measured evolution between Abell 1689 and 
Coma, decreasing A/3 by around 0.25 mag; removing the 
BCGs increases A/3 by around 0.05 mag. 

Fig.[6]shows the FJR measured for AbeU 1689 in F625W 
compared to that measured for Coma in the r'-band. Both 
de Vaucouleurs and Sersic surface photometry are shown for 
Abell 1689 but only de Vaucouleurs photometry is available 
for Coma. We have corrected the data points for size evo- 
lution and over plot the best fit DLM (omitting BCGs) in 
each case. The DLM fits the data well and both the Coma 
and Abell 1689 galaxies populate the same {m, ct} parameter 
space. 



3.3 Kormendy relation 

Table [7] presents the results for the fits to the KRs (Eq. [2| 
of Coma and Abell 1689 using LMs and the DLM, as was 
done for the FJR above. 
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Table 7. Parameters and marginalised uncertainties for the fits to the Faber-Jackson and Kormendy relations for Coma and Abell 1689 
as well as the CMR relation for Abell 1689. Model abbreviations are: MM=Mixture model; LM=Linear model; DLM=Double linear 
model. The column headed SEC shows whether a size evolution correction was applied f ij2.8.3l l. The Profile column shows if the surface 
photometry (COG) was calculated using de Vaucouleurs (dV) or Sersic (S) profiles. Note that in the double linear fits made to Coma 
and Abell 1689 data, the Coma data is archival and was measured using de Vaucouleurs COGs while for the Abell 1689 data we had 
the option to fit Sersic or do Vaucouleurs profiles. Where BCGs were removed from the fitting procedure, we removed galaxies #584 and 
#610 from the Abell 1689 sample and GMP2921 and GMP3329 from the Coma sample. For the KR and FJR, results in bold are the 
ones we consider to be most reliable: note that without size evolution corrections the FJR and KR luminosity evolutions disagree {2a). 
Taking the average of (i.e. combining the samples for) the KR and FJR results corrected for size evolution while excluding BCGs gives 
A/9 =0.22 ± 0.11 mag (see il3.5|l which we compare to stellar population models in i|4.2l 
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Figure 6. The de Vaucouleurs (top) and Sorsic (bottom) Kormendy (left) and Faber- Jackson (right) relations for Abell 1689. Squares 
represent good data while diamonds and arrows represent bad data (poor COG fits, or upp er limits on velocity dispersions) that were 
excl uded from the fits. Smal l black squares represent Coma data from lJorgensen et al.l l ll995al . in the r-band, corrected to the AB system) 
and |j0rgensen et al.l l ll999bl ) while large coloured symbols represent data for Abell 1689 (from HST/ACS F625W imaging and GMOS-N 
spectroscopy). Error ellipses on the Abell 1689 data are shown in light grey and account for correlations between each axis. The error 
bars are larger for the de Vaucouleurs data because they include the systematic errors associated with fixing the Sersic index to n = 4. 
Solid lines show the results of fitting the double linear models without BCGs: the light solid line represents the Coma data while the 
bold solid line represents the Abell 1689 data. Dashed lines illustrate the magnitude cuts (50% completeness) for the Abell 1689 sample 
(bold dashed) and the Coma sample (light dashed). Note that bad points (diamonds) may lie outside the plotting window. Galaxies for 
which o- < 50km s"! arc shown as upper limits in the FJR (#724, #884 & #906). 



The values for A/3kr in Table [T] like A/3f.th, were calcu- 
lated both at logi?e = 0.5 from individual LMs, and directly 
from the DLM. As before, we favour the values calculated 
by the DLM. As for the FJR, galaxies in Abell 1689 are 
fainter in the F625W band than Coma galaxies in the r'- 
band (A/3 > 0). Unlike the FJR, there is a significant differ- 
ence between the evolution measured from the LMs and the 
DLM: the KR is sensitive to magnitude cuts and selection 
effects which is unsurprising given that the cut runs almost 
parallel to the KR (unlike the FJR). Correcting for size evo- 
lution increases A/3kr only by around 0.05 mag, while in- 
cluding the BCGs has a negligible effect: clearly the KR is 
relatively robust to these two factors. 

Fig.[6]shows the KR measured for Abell 1689 in F625W 
compared to that measured for Coma in the r'-band. Both 
de Vaucouleurs and Sersic surface photometry are shown 
(for Abell 1689) and as before, we have corrected the data 
for size evolution and over plot the best fit DLM (omit- 



ting BCGs). The galaxies in Abell 1689 do not populate 
the same {(At)e, fie} parameter space as the Coma galaxies: 
there appear to be more bright, compact galaxies in Coma. 
Recall that we have corrected the Abell 1689 galaxies for 
size evolution, thus making them larger and fainter, but we 
did not apply any cuts in _Re or (/i)e a priori; the only cut 
was in absolute magnitude, which is well matched to the 
absolute magnitude cut in the Coma sample. Indeed, using 
SExTRACTOR effective radii as an estimate for size, there is 
no difference between the size distributions of our spectro- 
scopic sample and the parent sample with r' < 19.5. Thus 
the different locations of the Coma and Abell 1689 data 
along the KR are not due to size selection effects. 
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3.4 The effects of size evolution in the KR and 
FJR 

Without correction for size evolution and excluding BCGs, 
the measured A/3kr and A/?fjr disagree at the 2a level. 
Although the need to correct for size evolution is at a modest 
statistical significance, reviewing the effect reveals why the 
FJR was more affected than the KR. 

If we consider both the corrections to Re and (/i)e, us- 
ing our knowledge of how these two are related (Eq. \^ we 
can infer how they shift the whole relation up or down in 
{/i)e, which can be (mistakenly) interpreted as luminosity 
evolution of the stellar population. Using Eq. [T5]and Eq. [T7] 
in Eq. [2] we find 



^ie,z 



aKn\ogRe,z + /3kr + A/?K1 



(19) 



with 



A/?;,H = (aKR-5)Clog(l + z). 



(20) 

This A/3kr appears as a change in the surface brightness; it 
arises just from scaling the size of the galaxies and not from 
changes in the stellar populations. We see that the effects 
from changing Re and {/i)e counteract with each other in 
the KR: the typical value for qkr is ~ 4, which with C = 1- 
gives A/3kr = ^log(l + 2). This tells us that if we were to 
compare {/i)e vs. \ogRe for Coma to {fJ.)e,z vs. log_Re,z for 
Abell 1689, even without any changes in the stellar popu- 
lations, we would find that the Abell 1689 galaxies appear 
brighter than the standard KR prediction by log(l -|- z). 

Repeating this analysis for the FJR relation using Eq. 
I18l in Eq. [1] we find a similar relation to before 



M = «F.IRl0g(0-^) + /3f.ir + A/3f, 



(21) 



but this time with 



A/Jf.ir = -0.5aF.iR7?log(l + z). 



(22) 

As with the KR, one could easily misinterpret A/3fjr as 
luminosity evolution of the stellar populations, but it is in 
fact just an effect of size evolution. However, unlike the KR, 
the change is quite substantial. Our typical Ofjr is ~ —8 and 
with ri = 1, gives AM = 41og(l + z). If we were to compare 
M vs. logcr for Coma to Mz vs. logcr2 for Abell 1689 and 
there was no change to the stellar populations, we would still 
find Abell 1689 galaxies to be fainter than the standard FJR 
prediction by 41og(l + z). The canonical value for ofjr is 
~ 10 (1/ ~ CT*), which increases the effect to 51og(l -I- z). 

Thus size evolution has opposite effects in the KR and 
the FJR: it makes the KR of earlier galaxy populations ap- 
pear brighter, the FJR fainter, and we observe these changes 
in A/3fjr and A^kr. 



3.5 Luminosity evolution from the KR and FJR 

After correcting for size evolution, the differences between 
the KRs and FJRs of Coma and Abell 1689 (A/3kr and 
A/3fjr) agree. To avoid propagating multiple calculations of 
A/3 any further, we now consider which values are our best 
estimates. 

The values of A/3 calculated using the DLM account 
for magnitude cuts and selection effects and are fully 
marginalised over the other model parameters; the values 
calculated from separate LMs to Coma and Abell 1689 are 



only marginalised over the model parameters fitted to each 
cluster, depend on where (on the x-ax;is) we compare the two 
clusters, and do not account for magnitude cuts or selection 
effects. We favour calculating A/B using the DLM for these 
reasons. However, the DLM, as currently defined, cannot 
address downsizing. 

We should only consider values of A/3 which were cal- 
culated with size evolution corrections: without such correc- 
tions, A/3kr and A/3fjr disagree. Including BCGs reduces 
A/3fjr and has a negligible effect on A/3kr, but the changes 
are small and within the uncertainties. The average A/3 with 
BCGs is 0.23 mag while the average without BCGs is 0.22 
mag. We somewhat arbitrarily choose to use the values of 
A/3 calculated without the BCGs (see ^O)) . Finally, A/3 
measured using Sersic photometry in Abell 1689 is on aver- 
age slightly larger than that measured using de Vaucouleurs 
photometry for Abell 1689, but the difference is far less than 
the uncertainties. Consequently, we take the average of the 
KR and FJR offsets (excluding BCGs) using both Sersic and 
de Vaucouleurs photometry as the overall luminosity evolu- 
tion. We now seek the uncertainty in this value. The {fi)e 
and Re data used in the KR are strongly correlated to the 
magnitudes used in the FJR so we cannot treat the uncer- 
tainties in A/3kr and A/3fjr as independent. We combine the 
MCMC samples for A/3kr and A/3fjr from the fits to Sersic 
and de Vaucouleurs surface photometry, thereby assuming 
equal evidence for all these models and measure the differ- 
ence in luminosity between the Coma galaxies (measured in 
rest-frame r') and the Abell 1689 galaxies (measured in the 
observed-frame F625W) to be 0.22 ±0.11 mag. We do not 
attempt to quantitatively justify the assumption of equal 
evidence as the results from the different models are very 
similar. 



4 STELLAR POPULATION MODELS 

As discussed in i|T] changes in the scaling relations can be 
understood in terms of evolution of the stellar populations. 
In order to infer such evolution from the observations, it is 
necessary to make use of stellar population synthesis models 
which reproduce the SED given some assumptions regarding 
that population, such as the age, the IMF, the metallicity (Z) 
and the star formation history (SFH). The latter describes 
the distribution of star formation over time and in our case 
is assumed to be a delta function at some previous time 
(i.e. a simple stellar populations or SSPs), so has a single 
characteristic age. 

It is beyond the scope of this paper to perform an ex- 
haustive comparison of all the possible models now avail- 
able. We use two common p opulation synthe sis models 
(|Bruzual fc CharlotI |2003| and iMarastonI 120051 , hereafter 
BC03 and M05) to study the evolution of the scaling rela- 
tions. Both provide SSP SEDqj spanning a variety of ages, 
metallicities and IMFs, although we only consider models 
using a Salpeter IMF. 

We use these models together with the filter curves pro- 
vided by BC03 (SDSS and Gunn bands) and the GEMINI 



** http://www.cida.ve/~bruzual/bc200 3 (BC03) 

|http://www. leg. port. ac.uk/~inaraston/ (M05) 



and 
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Observatorjlfl (GMOS-N r' and g') and the Vega and Solar 
spectra provided by the HST calibration database (CAL- 
SPEQj). We make use of the COSMOLOPMJ python library 
to calculate look back times and redshifts given our cos- 
mological parameters (see Q. The model magnitudes are 
calculated so as to be comparable to observations that are 
K-corrected for bandpass effects {Kb), but not for colour 
effects {Kc)- 

Observations of objects at different redshifts equate to 
observations in different frames of reference (observer, rest), 
while all the model calculations are performed in the rest 
frame. This, together with the different magnitude systems 
(Vega, AB, Gunn) and the different epochs at which we ob- 
serve our galaxies, need particular attention. Consequently, 
where necessary, we label the specified filter with the mag- 
nitude system, the frame of reference and the epoch of the 
observation in parenthesis. For example, F625W(AB, z=0.0, 
T) refers to AB magnitudes of local objects in the F625W 
filter, while F625W(AB, z=0.183, T-2.26 Gyr) refers to ob- 
servations in F625W for an object at z=0.f83 (the look back 
time for an object at z=0.f83 is 2.26 Gyr for our adopted 
cosmology). The difference F625W(AB,z=0.f83,T-2. 26 Gyr) 
- r' (AB,z=0.0,T) is what we measure when comparing our 
Abell f689 data to the Coma data (i.e. A/3). We also find it 
useful to quote model values for which we do not have obser- 
vations: F625W(AB, z=0.183, T) refers to observations in 
the blueshifted F625W filter at the current age of the Uni- 
verse and so F625W(AB, z=0.183, T)-F625W(AB, z=0.183, 
T-2.26 Gyr) is just the luminosity evolution between now 
and 2.26 Gyr ago in the blueshifted F625W filter bandpass 
(approximately equivalent to the V-band). 

These stellar population models can be used in many 
ways, but we choose two distinct and very different tech- 
niques for the CMR analysis and KR/FJR analysis. We 
present the BC03 and M05 results for modelling the passive 
evolution of SSPs as a function of age and various metallic- 
ities in both cases. 



4.1 Analysis of the g'-r' CMR 

As discussed in ijl] the intrinsic scatter in the CMR tells us 
about the SFH of the galaxies. We can perform a similar 
analysis to BLE92 for Abell 1689, but using state of the art 
population synthesis models. Fig. [7]shows the rate of change 
of (g'-r') as a function of the age of the stellar population for 
the BC03 and M05 models with solar and super-solar metal- 
licities. Like BLE92 and BKT98, we smoothed the models; 
the resulting systematic error was no more than 0.1 mag. 
Our SSP tracks in Fig. [7] are broadly similar to those of 
BLE92 and BKT98, but the M05 models show a more rapid 
fall in the rate of change of colour before <4 Gyr, followed 
by very slow or even negligible further reddening beyond > 4 
Gyr. 

Following BLE92 and BKT98, we model the formation 
of galaxies with a uniform distribution of random events 
distributed over a time interval At with mean age (now) 



^ http://www.gcirLini.edu/sciops/instrumcnts/gmos/ 

^ http://www.stsci.edu/hst/observatory/cdbs/calspcc.htmI, 

post Feb. 2010 update. 

^ http://roban.github.coni/CosnioloPy/ 



if. We parameterise the time interval At as a fraction, /3, 
of the time available up to the end of At; when galaxies 



form between times t. 



U 



and 



At = /3tatop with < /3 < 1 (when /? = 1, galaxy formation 
is continuous from the start of the Universe to some final 
time tgtop so we have the largest possible scatter in ages and 
no synchronisation) . We can write At in terms of t/ , 



At: 



tuniv 



/3-1 - 1/2 



(23) 



where tuniv is the age of the Universe now. 

As noted by BLE92, if the slope of the CMR is produced 
by variation in metallicity with luminosity and all galaxies 
have have identical IMFs, then the RMS scatter of the CMR, 
o"cMR is related to the RMS scatter in the formation ages of 
the RS galaxies, (Jagc by 



CTcMR ^ 9(g' - r') 



dt 



= ./(age) 



(24) 



where the RMS scatter a.^^^ ~ At/3.5 for the uniform dis- 
tribution considered here (see BLE92 & BKT98). 

The observed scatter in the CMR can therefore tell us 
about the scatter in the formation times (synchronisation, 
P). But as galaxies become redder with time and the rate of 
change of colour decreases, a small scatter in the CMD could 
also imply a very old population with no synchronisation. 
Eq. [24] can be written 



a(g' - r') 
dt 



3.5crci 



p-^ - 1/2 



-ff 



(25) 



For a given /3, there is an age, tf, which is compatible with 
the observed scatter (Tcmr. fn Fig.[7l we plot the LHS of Eq. 
1251 (curves of 9(g' — r')/9t derived from stellar population 
models) and the RHS of Eq.[25]( "beta" curves) . Where these 
curves intersect gives the model age now that satisfies Eq. 
1251 While BLE92 assumed that the variation in the rate of 
change of colour with metallicity was sufficiently small that 
they could consider only one metallicity (in effect assuming 
that the slope of the CMR was zero), we investigate solar 
and super solar metallicities. Fig. [7] shows all these curves 
for the BC03 and the M05 models. Significant differences 
only exist for younger ages (i.e. for /3 ~ 0.6). We show the 
68% uncertainty limits for individual beta curves as dotted 
lines (based on the uncertainty in ctcmr). 

Fig. [7]shows that ctcmr in Abell 1689 is compatible with 
a large range of ages, depending on the degree of synchro- 
nisation. BLE92 and BKT98 found the same for Coma and 
that they argued that younger ages and smaller At are in- 
creasingly unlikely because such high SFRs would be easily 
visible but are not observed. However, we observe Abell 1689 
some 2.26 Gyr earlier than Coma: combining both results as 
though we had observed the same cluster twice over a 2.26 
Gyr interval adds another constraint: if the scatter in the 
CMR has not significantly changed over the last 2.26 Gyr, 
then the rate of change of colour also cannot have changed 
appreciably over 2.26 Gyr. To try and put an upper limit on 
how much the rate of change of colour could have changed 
over these 2.26 Gyr, let us consider a lower limit on the Coma 
CMR scatter of 0.036 mag (the scatter for just ellipticals 
within 600 Mpc in BTK98). The scatter of all RS galaxies 
2.26 Gyr earlier (in Abell 1689) is measured to be 0.054 mag; 
this gives a maximum rate of change of -0.008 mag/Gyr for 



Data and 2D scaling relations for galaxies in Abell 1689: a hint of size evolution at z^O.2 19 



0.25 




e 10 

Age NOW (Gyr) 

Figure 7. The rate of change of colour vs. the age of the stel- 
lar population for the gmos g'-r' colour using BC03 and M05 
models for solar and super-solar metallcities. The x-axis repre- 
sents tt+2.26 Gyr, which is the age now, rather than the age at 
z=0.183; this is to allow easy comparison with Fig. [8] 



the CMR scatter, or -0.028 niag/Gyr for the rate of change 
of colour. In Fig.[7l models that allow for such a small change 
over 2.26 Gyr have ages > 6 Gyr and P > 0.6. According to 
the exact intersections in Fig. [71 the galaxies in the RS of 
Coma and Abell 1689 formed between 0.55 < z < 1.55 (have 
5.5 < age (Gyr) < 9.5), if the Abell 1689 RS evolves into the 
Coma RS. The scatter could be overestimated in Abell 1689 
(we do not select galaxies by morphology or recession veloc- 
ity), so the upper limit on age, like the lower limit on the 
scatter, could be questioned; thus we are only confident that 
the galaxies in the RS of Coma and Abell 1689 are > 5.5 
Gyr old (formed at z > 0.55). Similarly, only models with 
P > 0.6 are consistent with this view. 



and not colour after 4 Gyr, according to panels (a) and (b). 
After 6 Gyr, the solar metallicity models show the slowest 
rate of change and there is significant difference between the 
BC03 and M05 solar metallicity predictions; because look- 
back time is not linearly related to redshift, the evolution 
with respect to redshift is almost asymptotically slow in (d) . 

In Fig.|8jc) and (d) we also over-plot the measured lumi- 
nosity evolution (A/?) found between the galaxies of Coma 
and Abell 1689 (see i33.5|l . The various stellar population 
models predict different ages for the Coma and Abell 1689 
galaxies from this measurement: considering all models and 
the la limits, the luminosity evolution between Coma and 
Abell 1689 is consistent with an age > 6.0 Gyr . Limiting 
the models to solar metallicities, the BC03 model suggests 
an age of 8.6±i;sj Gyr , while the M05 model suggests an 
age of 11.4±5'g'® Gyr (models only run to 15 Gyr): there 
is a discrepancy between the ages inferred from the BC03 
and M05 solar metallicity models, the latter being older by 
2.8 Gyr, although both lower limits are comparable. In the 
absence of any reasonable cause to discriminate between the 
BC03 and M05 models, we average the ages and uncertain- 
ties (i.e. combine the samples, assuming equal evidence for 
the different models) to conclude that the galaxies in Coma 
and Abell 1689 are now 10.2 ± 3.3 Gyr old (i.e. formed at 
z = 1.8±o°g). Note that because the M05 models are only 
calculated to 15 Gyr and the la upper limit on A/3 is greater 
than the model predictions at that age, we can only quote 
the upper limit from the BC03 model. 

It is customary to quote the luminosity evolution in the 
same band and our stellar population models allow us to 
calculate this. The BC03 and M05 solar metallicity mod- 
els both give a luminosity evolution of 0.25 mag in the 
F625W(AB,z=0.0) band and 0.23 mag in the r' (AB,z=0.0) 
band for their respective ages. 



4.2 Analysis of the KR and FJR 

Our principle goal is to determine what (if any) luminos- 
ity and colour evolution has taken place between now and 
z=0.183. We use the BC03 and M05 models to calculate 
the evolution of F625W(AB, z=0.183, T) - r(AB, z=0.0, 
T) as a function of the age, T (roughly a V-r colour). 
This is shown in Fig. EJa) and shows that after ~ 4 Gyr, 
F625W(AB, z=0.183, T)-r(AB, z=0.0, T) for the model pop- 
ulations changes slowly from 0.4 to 0.5 mag at 13 Gyr, so 
any ETC older than 6 Gyr now, would not be much bluer 
at z=0.183. Note that this colour is equivalent to the colour 
term of the K-correction in Eq. 1121 Kc'- there is variation 
with age and metallicity, which is why we do not try to 
apply it to our measurements. We calculate the luminosity 
evolution over 2.26 Gyr for F625W(AB, z=0.183, T) and 
r(AB, z=0.0, T) in Fig. [8jb): here we see the evolution is 
very similar for all models after 4 Gyr (to within ~0.1 mag) 
and drops rapidly with age. Finally, we calculate the lu- 
minosity evolution we measure when we compare the KRs 
and FJRs of Coma and Abell 1689 in Fig. M.c) and Fig. 
\Eld): the evolution of F625W(AB, z=0.183, 7-2.26 Gyr) - 
r(AB, z=0.0, T) is dominated by the change in luminosity 
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Figure 8. (a): the change in F625W(AB, z=0.183, T)-r(AB, z=0.0, T) colour as a function of age, T; (b): the change in luminosity for 
the r{AB, z=0.0, T) and F625W(AB, z=0.183, T) filters as a function of T; (c); the difference between F625W(AB, z=0.183, T-2.26 Gyr) 
and r(AB, z=0.0, T) as a function of T, line colours are as in (a), solid and dashed lines show the evolution and associated uncertainty 
inferred from the FJR and KR; (d): like (c), but as a function of formation rodshift, line colours are as in (a). All: magnitudes arc 
in the AB system; two different population synthesis models (BC03 and MOS) are shown for sub-solar, solar (Z=0.02) and super-solar 
metallicities. All models are consistent with galaxy ages > 6.0 Gyr . The BC03 Zq{=0.02) models suggest ages of 8.6±j| Gyr while 
the MOS Zq models suggest ages of 11.4±3 g' Gyr (models only run to IS Gyr). 



5 DISCUSSION 

We now discuss the context of our findings and highlight 
relevant caveats so that the reliability of our approach can 
be judged. We discuss the CMD first, followed by the Faber- 
Jackson and Kormendy relations and summarise the findings 
at the end. 



little or no synchronisation in their formation (/? > 0.6). We 
now ask if this result is in agreement with other clusters at 
other redshifts. Then we compare Coma and Abell 1689 to 
judge if we are making a fair comparison. We also discuss the 
usefulness of placing lower limits on the CMR scatter and 
finally explain why certain caveats in our approach should 
not strongly bias the outcome. 



5.1 The CMR of Abell 1689 

The scatter in the CMR of AbeU 1689, (Tcmr =0.054 ±0.004 
mag (for all galaxies with Mv < 17.9 mag within a projected 
radius of 570 kpc from the cluste r centre) is com p arable to 
the bi-weight scatter quoted by iTerlevich et al.l (|200ll ) of 
0.063±o;oi* mag. Combining these results, we break the /3- 
age degeneracy and show that the galaxies in both clusters 
have an age > 5.5 Gyr (and formed a,t z > 0.55) and had 



5.1.1 Abell 1689 in perspective 



L2l [Stanford et al.l (|l995h . lEllis et al 



As di s cussed — , _, ., . .. . 

|199^, [Stanford et al.l (|l998l ) and most recently iMei et al 
(t2009|) have studied the CMR scatter in clusters up to 
z ~ 1.3. The last two authors used a clipped bi-weight esti- 
mator to measure the CMR scatter; to our knowledge, we are 
the first to model and extract the CMR scatter using a mix- 
ture model which does not rely on arbitrary cuts, clipping 
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thresholds or parameter tuning. Furthermore, we present a 
CMR scatter that is marginalised over all reasonable slopes 
and intercepts, whereas previous work quotes values given 
a measured slope and intercept. That said, the findings of 
these papers agree with our work: the scatter in the CMR 
has remained relatively unchanged in the E TG population 
of massive gala xy clusters since z~ 1.3. Both [Stanford et al.l 
and lMei et al.l found the scatter in their clusters to be nearly 
always less than 0.1 mag with the average being around 0.06 
and 0.05 mag, respectively. Both authors were also careful 
to discuss how selection effects may bias the determination 
of the CMR scatter at higher redshift (so-called progenitor- 
bias): they selected galaxies morphologically to be Es and 
SOs, so mergers or spirals were excluded from the sample. 
We kno w that more (blue) spirals exist in clusters at higher 
redshift (jButcher fc Oemleiill984 ) so it is reasonable to con- 
sider that the population of galaxies on the CMR at higher 
redshift is not the same population observed at lower red- 
shift. Furthermore, clusters grow over time and more field 
ellipticals will join the cluster later on, producing a similar 
bias. But the bias is less here, because we do not morpholog- 
ically select our galaxies and at z=0.183 we are only looking 
2.26 Gyr into the past, which is less time for change in the 
ETG progenitors compared to Coma. 



5.1.2 Lower limits on the CMR scatter 

Until now, no one has tried to interpret a lower limit on the 
scatter of the CMR. We do so now tentatively, because we 
know that that scatter is likely to be smaller for just ellipti- 
cals and also for samples limited to more luminous galaxies. 
Taken at face value, the lower limit of the scatter measured 
here suggests that, with /3 — 1.0, the highest possible mean 
formation age for the RS galaxies of Abell 1689 is 9.8 Gyr in 
the past {z = 1.7, a la upper limit). This may be meaning- 
less because we know we have mixed populations (ellipticals, 
SOs etc.) which increase the scatter, but in future surveys it 
may be of interest to consider an upper bound on the mean 
formation redshift if we wish to dedicate effort to witness- 
ing it. If all the stars in the Coma ETGs formed in the first 
Gyr after the big bang (t, = 13.2 Gyr, /3 = 1.0), the CMR 
scatter would be ~0.005 mag - seven times smaller than 
that of Coma. It takes time for dark matter (DM) to cluster 
sufficiently to allow the baryons to start forming stars, but 
the difference between a scatter of 0.036 mag and 0.005 mag 
in the CMR of Coma implies there was a delay of ~3 Gyr 
between the big bang and the average SF period for rich 
clusters (assuming P — 1). This suggests that the majority 
of stars formed for z < 3 wh ich is in agreement with the 
cosmic star formation history (JMadau et al.lll996l . ll998h . 



5.1.3 Caveats 

BKT98 investigated exponentially declining SFHs truncated 
at different times and also looked at the effects of merging. 
They found the same conclusions as BKT92, validating the 
SSP approach used there (and here). We treat the random 
errors in magnitude and colour as independent, but in real- 
ity they will be correlated: noise in the central aperture of 
the r-band images will affect the colour term as well as the 
total r-band magnitudes. This will have minimal effect on 



our scatter measurement, but could slightly affect the un- 
certainties. We investigated the assumption that Abell 1689 
will evolve into Coma in H2.7.1I 



5.2 The Faber- Jackson and Kormendy relations 

After correcting for size evolution, magnitude cuts and selec- 
tion effects, analysis of the F JR and KR in Abell 1689 shows 
that the Abell 1689 galaxies in F625W are fainter than the 
Coma galaxies in rest frame r'-band by 0.22 ± 0.11 mag, 
which suggest and age of 10.2 ± 3.3 Gyr {z = 1.8±g°9). We 
now compare these findings with the literature and discuss 
the importance of correcting for size evolution. 



5.2.1 Evidence for the passive evolution of ETGs 

The luminosity evolution in the KR and FJR in Coma and 
Abell 1689 is consistent with a passively evolving stellar 
population from single burst 10 Gyr ago. This is in agree- 
ment with the long-standing findings that ETGs both lo- 
cally and at moderate redshift appear to be dominated by 
old passively evolving stellar populations that formed at 
z > 1. The majority of precision age determinations stem 
from measurements of absorption line indices: although the 
broad-band colours of the Coma ETGs have been known 
to be compatible with old populations for a long time 
(see Hl.l|l . the well documented age-met allicity degeneracy 
limits this approach for local clusters, but the degeneracy 
can be reduced by modelling the absorption line indices of 
the metals. Little work has bee n done specifically on the 
galaxy populations of Abell 1689: iMieske et all (|2004l ) stud- 
ies the ultra-compact dwarf fr action and lum i nosity func- 
tion of the cluster galaxies while lBanados et al.l (|2010|) studv 
the fa int end of the luminosity function. ICarretero et al.l 
(|2007|) studied line indices of ETGs in 4 massive clusters 
at z ~ 0.2 and found ages of around 10-15Gyr, but their 
sample did not include Abell 1689. Considerably more work 
has been done on the galaxy populations of Coma, in partic- 
ular with regard to determining the stellar p opulation ages 
JClcmcns ct al."2009'; Jorgcnscn ct al. 1999b: 'Mchle rtetal.1 
2003; Sanchcz-Blazqucz_ct al 2006, 2009; Rakos ct aQjioOTl ; 
Traeer et all l2008l: ISmith et all l2009l : [Harrison et"aLl |2010| . 



2011: Price et al. 12011 



There is still some uncertainty re- 
garding the ages (and presence of multiple populations), but 
the current evidence suggests ages between 5 and 10 Gyr for 
the Coma ETGs. Our age estimate from the broad band pho- 
tometry of Abell 1689 and Coma is therefore towards the top 
of this range, though the BC03 result alone (8.6±i;g Gyr ) 
is nearer the midpoint and the M05 result at the upper end 
(11.4±5"g'^ Gyr (models only run to 15 Gyr)); intriguingly, 
the age range implied by just the scatters of the CMRs in 
Abell 1689 and Coma (5.5 < age (Gyr) < 9.5) brackets the 
literature ages very well, even though earlier we expressed 
doubt regarding this upper limit (see M.iy 



5.3 The effects of size evolution 

The measured luminosity evolution between Abell 1689 

and Coma ag r ees with that of ot her clusters 

(iBarrientos et all Il996l: iPahre et al.l Il996l; IZiegler et all 
I2OOII : iLa Barbera et all l2003l : iFritz et al.l |2005| ). affhough 
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none of these studies corrected for size evolution. Most did 
not simulta neously study both the FR and FJR except for 
iFritz et alj who found that the galaxies appeared brighter 
in the KR than in the FJR (but the 0.07 mag difference is 
consid erably less than the 0.3 mag found herel. lSaglia et al.l 
lJ2010f ) allow for size evolution in their study of the FP 
out to z ~ 1. However, we are not probing a high redshift 
cluster, but a comparatively local cluster at z=0.183. The 
fact that we need size evolution to bring the KR and FJR 
into agreement suggests that whatever causes it is a gradual 
process and its effects are not limited to the high-z universe; 
indeed, the mechanism is probably still in action now and 
does not appear to have deviated from the high-z (1 + z)^''' 
laws. 

While there has been a reasonable passage of time be- 
tween z=0.183 and now (2.26 Gyr), observing the size evolu- 
tion taking place in the cluster environment and in the local 
Universe place important restrictions on models that wish 
to explain size evolution, which we now discuss. 

Merging: Wet or dry major merger s cenarios (e.g. 
Khochfar fc Silkll2006l : iHernguist eralll 19931 : iHopkins et d] 
20091 ) are unlikely to occur in the cluster environment be- 



cause cluster galaxies have a very high relative velocity 
(larger than the typical escape velocity of two massive galax- 
ies) making it prohibitively difficult to merge galaxies. An 
exception are central BCGs which would likely merge when 
multiple clusters merge as they lie at rest at the bottom of 
the potential. Minor o r late accretion of smaller galaxies (e.g . 
Naab fc Truiilldl2006l: [Mailer et al.ll2006l : iNaab et al.ll2009l : 



Hopkins et al.l |2009| . [2010|) could be more feasible because 



of the larger number of low-mass satellites, but in general 
the relative velocities would still be higher than the escape 
velocity unless the impact was at glancing angle. 

Accretion: Accretion of the intra-cluster stars (previously 
stripped from other galaxies) may be energetically plausi- 
ble, if no t in mass transfer rate s. However, whi l e the simu- 
lations of lBois et all (|201ll ') and lKhochfar et all (|201ll ') sug- 
gest t hat ATLAS3D slo w rotators grow from minor accre- 
tions, IStott et al.l (|201ll ) finds no evidence of size evolution 
in BCGs up to z < 1. 

Adiahatic e xpansion: The a ction of AGN expelling gas as 
proposed by iFan et al.l |2008| is perhaps unfeasible here be- 
cause we know that massive reservoirs of gas are stripped 
when galaxies fall into a cluster (observe d for gas rich spiral s 
falling into the cluster potential today. IChung et al.ll2009r ). 
However, the mass lost from evolving stars could be grad- 
ually stripped by the ICM. Adiabatic expansion could also 
occur when central mass (dark or luminous) is stripped or 
harassed to larger radii by high speed encounters, which are 
inevitable in clusters. 

Secular/Cosmological expansion: There are few secular ar- 
guments (the action of bars and/or resonances) for size evo- 
lution in the literature at present and perhaps these should 
be investigated further given our findings: if SOs are believed 
to be faded spirals (the vast majority of Es show fas t rotat - 
ing disc- like disk kinematics, lEmsellem et al.ll2007l . 1201 J ). 
then whatever action reduces disks to spheroids could also 
be responsible for the size evolution we observe. Recall also 
that size evolution in disks was predicted as a co nsequence 
of th em being truncated to the critical density (|Mo et al.l 
119981 ). which is entirely a secular/cosmological effect. 



5. 3. 1 Caveats 

Some caveats to consider for our KR and FJR analysis in- 
clude not correcting for colour gradients, invoking size evo- 
lution, not selecting galaxies morphologically or to be on the 
RS, and using non-independent PDFs as though they were 
independent. We have already investigated the assumption 
that Abell 1689 will evolve into Coma in ^2.7.11 

An important difference between the FJR and KR is 
that the latter is potentially sensitive to internal colour gra- 
dients; this may be a source of error in our work given that 
our F625W observations of Abell 1689 are approximately 
rest-frame V-band compared to the Coma r-band observa- 
tions. The colour gradients in ETGs cause them to be red- 
der towards the centre; this makes -Re larger at bluer wave - 
lengths. Using the Coma data from |j0rgensen et al.l (|l995al ). 
we find that _Re in Gunn g is a mean of ~ 5 ± 15% (or a me- 
dian of 2%) larger than Re in Gunn r. However, if one consid- 
ers the combined effects in the KR, we find that they cancel 
to first order. Using the same s.ps that gave rise to Eg. 1201 in- 
creasing Re by 5% leads to an apparent luminosity evolution 
of log(1.05) = 0.02 mag in the KR. If we calculate the differ- 
ence between the central col our ((At)e(g) — (/j ) e(r)) a nd the 
global colour (g-r) from the IJorgensen et al.l (|l995al ') data, 
we find a mean [(/i)e(g) — {M)e(r)] — [g — r] = of 0.05 ± 0.35 
mag (median of 0.07 mag). Allowing for the large scatter, 
this agrees with the above and is a negligible correction. 
Therefore we conclude that it is not necessary to correct for 
colour gradients. 

We have invoked the existence of size evolution to make 
sense of the different luminosity evolution seen in the KR 
and FJR. However, we have chosen a prescription for how 
the size and internal kinematics of the galaxies change with 
time based on the evidence to hand fi i2.8.3fl . Although this 
is the simplest approach, one could try to measure this pre- 
scription from the data by making the simple assumption 
that variation in the mass-to-light is determined only by 
passive evolution of the stellar population, is independent 
of galaxy luminosity (or stellar/dynamical mass) and is the 
same for the KR and the FJR. 

Readers may have noticed that we do not select 
Abell 1689 galaxies in the KR and FJR to be morphological 
E/SOs or to be on the RS (even though we have that infor- 
mation from the CMR analysis). This is deliberate: we do 
not want to bias ourselves away from including (blue, spiral) 
progenitors of ETGs in Coma. However, blue spirals tend to 
exclude themselves from our sample because they have ve- 
locity dispersions below our spectral resolution (most likely 
because the luminosity is dominated by light from central 
dynamically cold star forming regions). Of the 38 galaxies 
that have sufficient quality spectra and HST imaging to en- 
ter the fitting procedures, only one is found not to be on 
the RS (#435). Thus although we tried not to select just 
RS galaxies, we end up being significantly biased towards 
them. 

When calculating the difference between the evolution 
seen in the KRs and the FJRs, we treated the PDFs as 
though they were independent normal distributions. Neither 
of these assumptions is true (the Coma and Abell 1689 data 
used in both are obviously correlated), but it is the best 
approach we are reasonably able to take. 
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6 CONCLUSIONS 

We have presented the Kormendy, Faber-Jackson and 
colour-magnitude relations for ETGs in Abell 1689 using 
HST/ACS imaging, GEMINI/GMOS imaging and GEM- 
INI/GMOS spectroscopy and conclude; 

(i) The intrinsic scatter in the colour-magnitude rela- 
tion of Abell 1689 places degenerate constraints on /3 (the 
ratio of assembly timescale and time available) and the 
age of the population: specifically, it is consistent with 
the galaxies in the colour-magnitude relation having ei- 
ther formed randomly at high redshift [z >2), or at lower 
redshift with increasing synchronisation (smaller assembly 
timescales) . However, assuming the intrinsic scatter of Coma 
and Abell 1689 is that of the same cluster observed twice 
over an interval of 2.26 Gyr breaks this degeneracy and lim- 
its /3 to be > 0.6 (little or no synchronisation) and the age 
of the red sequence to be > 5.5 Gyr (formed at 2; > 0.55). 

(ii) After accounting for size evolution effects, the F625W 
Kormendy and Faber-Jackson relations both show a simi- 
lar change in luminosity compared to the rest frame r'-band 
relations of Coma: 0.22 ± 0.11 mag. This is consistent with 
passive evolution of the stellar populations from a single 
burst of star formation long ago: the galaxies in Coma and 
Abell 1689 have an SSP age of 10.2 ± 3.3 Gyr (i.e. formed 
at z = 1.8±o°9), which agrees with the CMR analysis above. 
However, not accounting for size evolution causes the Kor- 
mendy and Faber-Jackson relations to be inconsistent and 
they then disagree about the amount of luminosity evolution 
at the 2(T level. 

(iii) We therefore find weak evidence that size evolution 
appears to have taken place in the cluster environment in 
the last 2.26 Gyr; if true, this places interesting constraints 
on the models, favouring harassment or secular mechanisms 
over major or minor merger scenarios. 



APPENDIX A: THE CURVE OF GROWTH 
(COG) TECHNIQUE AND ASSOCIATED 
ERRORS 

The COG method we use is in essence very simple: locate 
the centre of the galaxy via some means (in our case, via 
ellipse fitting); integrate out in radius to create a curve-of- 
growth function; fit an analytical curve-of-growth (Eq. llOf) . 
varying the model parameters (apparent magnitude, m; ef- 
fective radius. Re and Sersic index n) to minimise the square 
of the residuals (x^)- 
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Figure Al. An azimuthal average of multiple HST/ACS PSFs 
(black points); radially binned points (red) and MGE approxi- 
mation (green). This MGE approximation is used with Ea. lA2l to 
account for PSF effects when fitting COGs. 



(|l993h . lTruiilloet al.1 |200li) and lTruiillo et all (|2001bh de- 
scribe methods to correct for this. However, the HST/ACS 
PSF is not well approximated by a single Gaussian or Mof- 
fat function, for which these methods were developed. We 
approx imate the AGS P SF with a multi-Gaussian expansion 
CMGE. lBendinellilll99ll '). We azimuthally average (4-fold re- 
flect) multiple PSFs from across the FOV (see ^2^211 . bin 
radially, and then use a ID MGE fitting routine to calculate 
a radial approx imation (made available by M. Cappellari 
ICappellarill2002l ). Fig. lAll shows an example. 

We make use of the well known result that the con- 
volution of the Sersic function I{R,Re,n) with a Gaussian 
G{r,a) is r{R,Re,n,a) = 



— - / I{x,Re,n)exp 



x^ +R' 



B(^] (A2) 



X 



E [^' - ^'^»]' 



(Al) 



where j? is a modified Bessel function of the first kind, zeroth 
order (JMoffatl [l969l : iBendineUi et al.lll982l 'l. Summing this 
according to the amplitudes of the MGE then provides a 
good approximation to the full PSF convolution. We do not 
account for ellipticity in the fitting process. A Levenberg- 
Marquardt algorithm (MPFItH) minimised x^ to fit Re, le 
and n (we fix n = 4 for a de Vaucouleurs profile). 

We find good agreement between this ID approach and 
Sersic profiles convolved with the full 2D HST/ACS PSFs: 
systematic errors in recovered parameters {Re,Ie, m, Rele'^) 
are always <1% over a wide range of parameter space. 



where the sum is performed over samples i for the observed 
COG Di and model COG Mi. We sample the COG at the 
plate scale of the input image. 

We now discuss how we deal with: PSF effects, contami- 
nation and masking of other nearby galaxies, and estimating 
our uncertainties. 

Al PSF effects 

The instrumental PSF of the image creates a systematic 
difference between the observed and true CQG. ISaglia et al.l 



A2 Contamination and Masking 

The galaxies in Abell 1689 are densely packed and many 
overlap in projection. To remove this contamination, we ini- 
tially bin azimuthally averaged profiles of the chosen galaxy 
and sigma-clip from the mode of each bin. If necessary, we 
further manually mask the images, ensuring identical masks 



http://www.physics.wisc.edu/~craigm/ 
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A3 Estimating errors (random and systematic) 

We use Monte-Carlo simulations to estimate uncertainties 
in the COG parameters. Generating Sersic profiles of vari- 
ous Re, le and Sersic index n convolved with a typical 2D 
HST/ACS PSF (see ig^J, we add random sky ofltsets and 
read and shot noise to simulate real observations. We use 
100 realisations for each parameter combination [Re , m and 
n) over the parameter space of the obervations. We also 
vary the radius at which we stop the COC fitting according 
to the statistics of the real data. This approach allows us 
to estimate random and systematic uncertainties (the latter 
always being < 1%). 

The simulations show that log_Re and log/e are very 
strongly correlated. Decomposing the covariance matrix into 
principle components, the first component always accounts 
for >98% of the variation: i?e-^e'^ is best constrained 

while _R~°'' /e is most poorly constrained. Knowing the 
covariance matrix as a function of input Re and m allows 
us to use it when fitting the scaling relations; treating Re 
and le as independent significantly overestimates the error 
in the direction perpendicular to the KR which is important 
when minimising perpendicular residuals (see ^Cjl. 



Figure A2. Examples of the masking procedure before fitting 
a curve of growth model, (a) HST/ACS F625W image of galaxy 
#286 which is relatively uncontaminated; (b) shows the mask ap- 
plied to this image (black pixels) and (c) shows the result of the 
image reconstruction from replacing masked pixels with a 2-fold 
reflection, (d) HST/ACS F625W image of galaxy #435 which 
has considerable contamination; (e) shows the mask for this im- 
age and (f ) shows the reconstructed image, produced by replac- 
ing pixels with a 2-fold reflection or azimuthally averaged values. 
(g) galaxy #635 is so badly contaminated that a reliable mea- 
sure of the COG is almost impossible; (h) shows the mask (leav- 
ing enough information to provide an azimuthal average) and (1) 
shows our best reconstruction. Although the reconstruction in (i) 
and the associated COG appear convincing, such results always 
score a 3rd class rating because of the large uncertainties. Scales 
on the X- and y-axis are in arcseconds, centred on the galaxy. The 
greyscale is logarithmic and flxed throughout, to highlight faint 
structure. 



for both the de Vaucouleurs and the Sersic COC fits to pre- 
vent systematic differences. To replace masked pixels, a sim- 
ple 2-fold reflection about the galaxy centre was usually suf- 
ficient. When masking was more severe, after a 2-fold refiec- 
tion, we azimuthally binned the image, averaged the modes 
of bins at each radius, and replaced masked pixels with this 
average. This reconstructed strongly contaminated images 
and gave good COG fits. However, as a method of quality 
control when the galaxy image was minimally contaminated 
and the model COG fitted the measured COG to within 2%, 
we rated the fit 1st class; when the contamination was again 
relatively minor but the model COC didn't fit the measured 
COC to within 2%, the fit was rated 2nd class; finally, when 
contamination was severe, such that we had little confidence 
in the measured COC and thus the model fit, we rated it 3rd 
class (shown in Table lU . Examples of masking and image 
reconstruction are shown in Fig. IA2I 



APPENDIX B: MATCHING SPECTRAL 
RESOLUTIONS 

The CMOS arc and sky lines appeared top-hat like because 
the slit width, not the grating, dictated the spectral resolu- 
tion. Matching spectral resolutions in this case is difficult: 
firstly, the size or shape of the objects being observed in 
the CMOS slits can define the resolution and spectral PSF 
if they are smaller than the slit width; secondly, the ker- 
nel required to match the spectral PSF of the stellar library 
with the CMOS spectra will not be Gaussian. Fortunately, 
the seeing FWHM (estimated by Gaussian fits to the pro- 
files of guiding/ alignment stars) is greater than the 0'.'75 
slit width (varying between 0'.'74 and l'.'49 with an aver- 
age of l'.'06). Thus the illumination of the slits was roughly 
uniform and similar to that of the arc or sky lines, which 
enables us to use arc or sky lines as a reference for the spec- 
tral PSF of the galaxy spectra. An additional complication 
in matching spectral resolutions is that the width of the 
CMOS arc lines decreases over the wavelength range used 
to extract kinematics; if one matches local stellar templates 
to the observer-frame arc spectra then when later redshift- 
ing the templates to the galaxy rest-frame, the non-uniform 
spectral resolution would introduce a systematic difference 
between the stellar and galaxy resolutions. 

We developed a technique to account for all these ef- 
fects by finding a transfer function which, when convolved 
with the spectral profile of the stellar library, reproduces 
the spectral profile of the galaxy spectra (which is equiv- 
alent to reproducing line profiles of a blueshifted arc spec- 
trum). We found the CMOS arc lines were w ell described by 
a truncated Causs-Hermite ( CH) expansion iJGerhardll 19931 : 
Ivan der Marel fc Fran^ll993l '): 



T(A) = 



7 



2Tva 



exp 



(A - Aq) 

2cr 



J^h^H^H, (Bl) 



which is a Gaussian centred on Ao with scale factor 7, disper- 
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sion a weighted by a sum of Hermite polynomials Hi. Taking 
/io = 1, we quantify deviations from a Gaussian using the 
truncated series {/!,2, /14, /le} (odd Gauss-Hermite moments 
are not required to describe a symmetric shape). Further- 
more, given that positive values of hi (where i is even) pro- 
vide more top-hat like profiles, whereas negative values give 
profiles with large wings and sharp peaks, only positive hi 
values are required. 

We then saught a Gauss-Hermite expansion which, 
when convolved with another Gaussian of FWHM equal 
to that of the stellar library resolution (lA), matched the 
blueshifted arc lines. We accomplished this by minimising 
the x^ difference given by the expression 



J2[A',-{S.{a')^T,{a,h,))] 



(B2) 



where A' is the (blueshifted) arc line being fit, 5* is the 
spectral profile of the stellar library (width a'), T is the 
Gauss-Hermite transfer function (Eq. IB1|) and the sum is 
performed over pixels. Fig. |B2] illustrates the parametric re- 
sults of this process at various wavelengths (and what would 
happen if we were not to account for the redshift effect). It 
is evident that a single transfer function with parameters 
{a,h2,hi,h6} = {0.4281,0.676,0.0,0.435} is sufHcient to 
degrade the spectral resolution of the stellar library to the 
blueshifted GMOS arc spectrum. We did occasionally expe- 
rience a problem of multiple minima when fitting the con- 
volved Gauss-Hermite transfer functions but careful choice 
of starting parameters solved the problem. Note that be- 
cause the width (A) and properties of the transfer function 
are uniform in A space, we convolved the stellar library with 
T before resampling it to log — A space. Fig. IB2I compares 
the profile of the matched stellar library to a typical profile 
for a (blueshifted) GMOS arc line: they are barely distin- 
guishable from each other, confirming we are able to match 
the spectral profile of the stellar library to that of the galaxy 
spectra. 
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Figure Bl. The best-fit transfer function, matching the spec- 
tral response of the stellar library to GMOS arc emission lines at 
various wavelengths. Parameters for the best-fit transfer func- 
tions (Eq. IBll l to blueshifted arc spectra are approximately 
uniform for all wavelengths (asterisks) and are well approxi- 
mated by a single function with parameters {cr, /12 , /i4 , /le } = 
{0.428A,0.676, 0.0, 0.435}. The widths of the GMOS arc lines 
(cr, black diamonds) are not constant (in A or log — A space), 
which can introduce a systematic error when stellar templates 
are redshifted to the recession velocity of the cluster. The cor- 
rect transfer function parameters, accounting for this effect are 
shown as green diamonds (previously the black asterisks fitted 
to blueshifted arc spectra, but now redshifted back). The widths 
of the transfer functions if we had not corrected for the redshift 
effect are shown as red diamonds: these appear roughly equal to 
the raw arc widths (black diamonds) but this is bad because the 
transfer function still needs to be convolved with the response 
of the stellar library; after this, the stellar templates would have 
a larger spectral resolution than the galaxy spectra, introducing 
a systematic error into the kinematics (dispersions would be too 
low). 
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Figure B2. A fiducial GMOS arc line (black) with the raw stellar 
library instrument profile (blue), best-fit Gauss- Hermite transfer 
function (red) and matched stellar library profile (green) over 
plotted. Note that the matched stellar library profile is virtually 
identical to the GMOS arc line profile. 
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while the combined projected variance perpendicular to the 
linear model is 

(C7) 



(C8) 



Equally, we could have introduced the rotated coor- 
dinate system of the linear model, Z' = [a::',?/'] which de- 
scribe coordinates parallel and perpendicular to the regres- 
sion line. They are related to the non-rotated coordinates 
by [a;', y'] = R. [x, y] — /?' where R is the usual rotation ma- 
trix which depends only on 9. There is an obvious equality 
between Ai and y-. There is no need to persist with this 
rotated coordinate system here, but it will be useful in the 
next section. 

The likelihood of observing all our data points given a 
particular choice of linear model is then the product of the 
frequency distributions, 



= 11^(2^* 



Xi , (Jx 



,<^v,i,o:,l3,ai). 



(C9) 



APPENDIX C: FITTING SCALING 
RELATIONS USING AN MCMC APPROACH 



In i|2.9l we briefly discussed the different models we use to 
fit the scaling relations in Coma and Abell 1689. Here we 
discuss those models in mor e detail. We have ta ken on board 
the recent remarks made bv lHogg et al.l (|2010l ) and we turn 
to Markov Chain Monte Carlo (MCMC) methods to 'fit' the 
data. We use three different models; each one is discussed 
below. 



CI The (Simple) Linear model 

This model is described by iHogg et al.l (|2010l ) in their Eq. 
35. We assume that our data {xi,yi} with (normally dis- 
tributed) uncertainties {(Txi,o-yi} and covariance a^yi — 
pxyiO'xiCTyi come from a linear relation y = ax + /3, with 
intrinsic scatter ai. Given a data point Xi with uncertainty 
axi we can define the expected frequencjij of observing the 
corresponding data point in the range [j/i, yt + y] with uncer- 
tainty ayi given some particular model parameters {a, /3, cri} 
as 



p{yi\xi,crxi,ayi,a,l3,ai) = 



1 



VM< + ^'i 



: CXp 



where 

A, = 



/?' 



sK + 'x?) 

(Cl) 
(C2) 



measures the perpendicular distance between the linear 
model and a data point because 

V — [cos 6, sin S] 
/3' = /3 cos 9 
sin 9 



cos 9 



(C3) 
(C4) 
(C5) 

(C6) 



Standard fitting procedures based on minimising x^ or 'max- 
imising the likelihood' then vary the model parameters to 
maximise this expression and find the mode of the likeli- 
hood distribution. However, the Bayesian view on this is to 
say that the posterior probability distribution of the model 



p(a,/3,o-i|{y,}i=i,Q) 



C.pia,m) 
P({j/J£ilO) 



(CIO) 



see lHogg et~aLl 1I2OIOI '). note #13 



where Q is shorthand for all the other (not always prior) 
knowledge of the problem, such as the Xi, {f7x,i,(yy,i} etc. 
Readers familiar with Bayes Theorem will recognise all the 
terms above, but for the unacquainted, p{a, f5\Q) contains 
additional information about the problem (and is referred to 
as the prior) and p{{yi}iLi\Q) is the normalising constant 
we evaluate by insisting that our final posterior distribution 
is normalised to one. Note that maximising the likelihood is 
equivalent to maximising the posterior given uniform priors 
on all parameters. 

MCMC samplers can probe the above likelihood with 
different model parameters and provide a set of sample 
points that are distributed according to the posterior prob- 
ability distribution (PDF) because they were drawn from it. 
This is our preferred method here, for two reasons. Firstly, 
it relaxes the standard assumption that the parameters are 
normally distributed; while normality may be a good ap- 
proximation for m and b, it is certainly not so for or be- 
cause f7i > 0, providing a hard (asymmetric) edge to the 
allowed distribution which we can set as a prior. Secondly, 
using MCMC here in a relatively simple model prepares us 
for more complicated models. 

In terms of im plementation, we use the PYMC library 
(jPatil et al.ll2010l ) to sample parameters 9,l3',ai using the 
standard Metropolis-Hastings MCMC algorithm. We choose 
relatively unrestrictive priors that are flat but always over 
a finite range: we let 9 range between 9o ± 7r/2 (where 9o is 
a good initial guess; this avoids searching for solutions near 
the boundaries); the intercept (/3') is limited to a finite range 
because we can't sample to ±cxd; the intrinsic scatter (ai) is 
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not allowed to be negative and also has some sensible upper 
limit. 

Whenever we quote a best-fit or most probable value for 
one of the parameters, we quote the median of the result- 
ing samples (which by default is marginalised over all other 
parameters). The MCMC usually requires reasonably good 
initial guesses to converge rapidly and such a guess can be 
provided by the 'maxi mum likelihood' e stimates or by eye. 

As highlighted bv lHogg et al.l (|2010l ). outliers can cause 
havoc in this simple linear model and iJC3] discusses our im- 
plementation of the solution proposed by those authors. 



C2 The (Complex) Double Linear Model 

It is often the case when studying scaling relations that one 
wishes to compare the intercept or offset of one relation to 
another. A common solution is to calculate the offset at a 
particular value along the x-axis, but such a result is depen- 
dent on the different slopes assumed for the two relations 
(as well as their intercepts). Alternatively, if there is good 
reason to believe the slope of the relations hasn't changed, 
one might fit the slope of the local relation (which usually 
has more data) an d then impose thi s fixed slope on the high 
redshift data (e.g. iBarr et al.ll2006l ): then at least there is 
no error created by imposing two different slopes, although 
there is still a strong dependence on the slope itself, the 
intercepts and any internal scatter being derived. 

To circumvent these problems, we propose a model con- 
sisting data drawn from two populations, Zj^ — [a:^i,i, yi,i] 
and Z^i — [2:2,1,2/2,1], each of which is described by a lin- 
ear relation (as in i ]Cl|l but we constrain the linear relations 
to have the same slope a and internal scatter oi; we then 
let the first population have intercept /3i while the second 
population with intercept /32 is offset to this by A/3, such 
that (32 = /3i + A/3. In this way, the two linear models share 
the same properties, except for the intercept and we have 
parameterised the variable we are most interested in as A/3. 
The likelihood for this model thus far would be 



i.'DL 



i i 



where 
Pi.i = 

P2,j = 



/ 



27ra?_, 
1 



: exp 



27ra2^^. 
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(Cll) 

(C12) 
(C13) 



give the frequency distributions of the individual data points 
for the two populations: Ai,i=v .Zi,i — /J( defines the resid- 
uals for the first population and A2,j similarly defines the 
residuals for the second population; the expressions 



2 2 

= cri,di + <^l 

2 , 2 

= cr2,dj + O"! 



(C14) 
(C15) 



give the total variances perpendicular to the relation. Hope- 
fully the definitions of gi, (Ji di and 02,di are obvious to the 
reader from iJCll 

One would also like to account for magnitude cuts and 
selection effe cts which are k n own to cause bias i n the derived 
parameters (|Willickl 1 1994 iTeerikorpil 1 19971 : ISagha et all 



I2OOII : iLa Barbera et al.ll2003l : iButkevich et al.ll2005^ . In the 
KR and FJR, magnitude cuts have the effect of masking a 
region of parameter space from being studied. In both cases, 
the division between the unsampled and sampled regions can 
be defined by a linear relation of the form y — a^^tX + Pant- 
For the KR, 

(C16) 
(C17) 



Qcut = 


= 5 




/3cut = 


= m<,„t + 2.51og(27r) 






-7.51og(l + z)- 51og 


^ kpc 
i arcsec 



while for the FJR, 

Qcut = (C18) 

/3cut = m + 2.51og(l + 2)-51ogDL + 5 (C19) 

The ratio of the spectroscopic samples to the parent 
samples as a function of m for both Abell 1 689 and Coma 
can b e fit by a selection function of the form (jWegner et al.l 
Il996h 



^'-2 



erf 



nicu 



V^5„ 



(C20) 



where erf is the error function, m^^t is the position of 
the cut (also the 50% completeness limit) and 5,n is the 
width/uncertainty of the cut (recall that Eg. IC20l describes 
a Heaviside step function at rricut convolved with a Gaus- 
sian with a — Sm)- We fit the selection function (Eq. IC20|| 
to the spectroscopic samples of Coma and AbeU 168 9 fol- 
lowing the approach described in lWegner et al.l ([199a). For 
Abell 1689, using SExtractor magnitudes, we find that 
the spectroscopic sample is consistent with selecting 60% 
per magnitude bin of the parent sample down to nicut = 19.5 
mag with uncertaint y 5m = 0.62 mag; for Co ma, we find the 
spectroscopic data of ljorgensen et al.l (119961 1 consistent with 
sampling 100% of the population down to rricut = 14.8 (AB) 
mag with uncertainty Sm = 0.5 mag. 

We model t he magnitud e cuts of t he two samples in 
a similar way to ISaglia et al.l (|200ir i and iLa Barbera et al.l 
(2003! ) by scaling the frequency distributions (i.e. the pi) by 
a factor 1/fi to account for the parameter space masked by 
the cut. This requires us to make an assumption about the 
distribution of the data along the regression line which we 
assume is uniform, starting at x'^^^ — inm{x'i) and ending at 
2;max = max(a::';) thus having width L = max(a:';) — min(a;i) 
(which is not a free parameter, but is set by the data and 
values of a and /3) . Therefore the probability space masked 
by the magnitude limit as a fraction of the total probability 
space is 
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(C21) 



where a^ is again the total variance perpendicular to the 
relation and we have evaluated the integral in the rotated 
coordinates of the linear model, {x',y'}. Eq. IC21I can be 
solved using the integral of the error function and substitut- 
ing the relation y' = a'^^^x' -\- 13'^^^ where {Ocut, /3cut} describe 
the magnitude cut in the rotated coordinate system. Know- 
ing ^cut = arctan(acut), it follows that 

aLt = tan(6',„t - 0) (C22) 

^cut = /3cut [cos e^^t ~ tan(6l,„t - 6) sin 6><,„t] - P cos t(C23) 
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Thus Eq. [C2T1 becomes 



/. 



+ 



PL 



erf 
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-— exp 
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(C24) 



(C25) 



Unlike ISaglia et alj and lLa Barbera et all we do not assign 
zero probability to objects with m > m^ut, but we scale pi by 
the selection probability Si. This accounts for the gradual 
nature of our magnitude cut. 

For the Coma sample, there is no other obvious selec- 
tion effect and Si accurately describes the ratio of the spec- 
troscopic sample to its parent distribution as a function of 
magnitude. However, the Abell 1689 data does suffer fur- 
ther selection effects which are not related to the magnitude 
cut above and should be corrected for. Firstly, only 60% of 
galaxies were selected per magnitude bin on average. Fur- 
thermore, we sample more than 60% of the most luminous 
galaxies due to low number statistics. Finally, we have the 
constraint that the data which enters the KR and FJR must 
have HST surface photometry, which further limits galaxies 
to be in the core of the cluster. Because of mass segregation, 
this biases the sample to more massive (luminous) galax- 
ies. We calculate the ratio of the ideal sampling (given by 
Si) to the actual sampling as a function of m which can be 
used as a w e ight Wi to correct for sampling effects, as in 
ISaglia et al.l (|200l|). Thus, accounting for the cut in mag- 
nitude with Si and fi, and additional selection effects as 
a function of magnitude with weights Wi, we can describe 
the frequency distribution of a single data point, given the 
model parameters to be 



Pi,i = 



SlA 



fi,iV2nai^i 



exp 






The final likelihood for the model is thus 



•LdlM — / , Pl,i + / ^ 



P2,: 



(C26) 



(C27) 



where the two sums sum over the data of the two different 
populations. 

Using the same MCMC set-up as before for this dou- 
ble linear model allows us to measure the marginalised off- 
set between two scaling relations while fully accounting for 
magnitude cuts and selection effects. 



C3 The Mixture Model 

We can also describe a scaling relation using a dou ble- 
Gaussian mixture model (Exercise 14 in iHogg et ali ) in 
which we again have data {xi,yi} with uncertainties 
{cxijO-yi} that come either from a linear relation {y — 
ax+P, model A, with intrinsic scatter ai) or another outlier 
distribution (model B with expectation value Yi, and intrin- 
sic scatter erf,) that overlaps that of the linear relation; thus 
bad data is modelled as coming from model B, whereas good 
data is modelled as coming from model A. The goal is to de- 
termine which data is good, and which bad. Say we have A'^ 
data points and we don't know a priori which distribution (s) 
the data points are drawn from, but we assume it's either 
A or B. Model A is described by a normal distribution just 



like the linear model (gCT 



PmodelA 



^2n{ai, + ai 



: exp 



A? 



2{al + ai) 



whereas model B is a simpler normal distribution 



PmodclS 



^2TT{al + a 



: exp 



(y« - Y,f 



2Ki+'-? 



(C28) 



(C29) 



Model B is described by just Yt and at, rather than any 
linear relation: it's just a scattering of points across all x, 
centred on Yb. 

We want to be ab le to reject data points when con- 
structing the likelihood. iHogg et al.l (|201G| ) explain how this 
is achieved by a summation over a mixture of the linear 
distribution (model A) and the outlier distribution (model 
BV3 such that the likelihood for the model is 



||[(1 — P(,).PmodclA +P6-Pmodcls] 
1 - Pb f A? 



(C30) 
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where pt is (a single parameter for) the probability of any 
given data point being bad. We again 'fit' the parameters 
9, P' ,ai,pt,,Yt,o-t, via MCMC sampling. The main advan- 
tage being a less subjective approach to removing outliers: 
we let the data and our knowledge of its uncertainty decide 
what to reject and what to keep. Again, the MCMC requires 
reasonably accurate initial guesses for the parameters to con- 
verge rapidly onto an acceptable result and such a guess can 
be provided by the 'maximum likelihood' estimates or by 
just being sensible. 
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